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Abstract. We introduce the classified stable matching problem, a problem motivated by academic 
hiring. Suppose that a number of institutes are hiring faculty members from a pool of applicants. Both 
institutes and applicants have preferences over the other side. An institute classifies the applicants based 
on their research areas (or any other criterion), and, for each class, it sets a lower bound and an upper 
bound on the number of applicants it would hire in that class. The objective is to find a stable matching 
from which no group of participants has reason to deviate. Moreover, the matching should respect the 
upper/lower bounds of the classes. 

In the first part of the paper, we study classified stable matching problems whose classifications belong 
to a fixed set of "order types." We show that if the set consists entirely of downward forests, there is a 
polynomial-time algorithm; otherwise, it is NP-complete to decide the existence of a stable matching. 
In the second part, we investigate the problem using a polyhedral approach. Suppose that all classifications 
are laminar families and there is no lower bound. We propose a set of linear inequalities to describe stable 
matching polytope and prove that it is integral. This integrality allows us to find various optimal stable 
matchings using Ellipsoid algorithm. By studying the geometric structure of fractional stable matchings, 
we are able to generalize a theorem of Teo and Sethuraman in the context of classified stable matching: 
given any number of stable matchings, if every applicant is assigned to his median choice among all stable 
matchings, the outcome is still a stable matching. Finally, a ramification of our result is the description of 
the stable matching polytope for the many-to-many (unclassified) stable matching problem. This answers 
an open question posed by Sethuraman, Teo and Qian. 
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1 Introduction 



Imagine that a number of institutes are recruiting faculty members from a pool of applicants. Both 
sides have their preferences. It would be ideal if there is a matching from which no applicant and 
institute have reason to deviate. If an applicant prefers another institute to the one he is assigned to 
(or maybe he is unassigned) and this institute also prefers him to any one of its assigned applicants, 
then this institute-applicant pair is a blocking pair. A matching is stable if there is no blocking pair. 

The above scenario is the well-studied hospitals/residents problem |7|llj in a different guise. 
It is known that stable matchings always exist and can be found efficiently by the Gale-Shapley 
algorithm. However, real world situations can be more complicated. An institute may have its own 
hiring policy and may find certain sets of applicants together unacceptable. For example, an institute 
may have reasons to avoid hiring too many applicants graduated from the same school; or it may 
want to diversify its faculty so that it can have researchers in many different fields. 

This concern motivates us to consider the following problem. An institute, besides giving its 
preference among the applicants, also classifies them based on their expertise (or some other criterion). 
For each class, it sets an upper bound and a lower bound on the number of applicants it would hire. 
Each institute defines its own classes and classifies the applicants in its own way (and the classes 
need not be disjoint). We consider this flexibility a desirable feature, as there are some research fields 
whose boundaries are blurred; moreover, some versatile researchers may be hard to categorize. 

We call the above problem classified stable matching. Even though motivated by academic 
hiring, it comes up any time objects on one side of the matching have multiple partners that may be 
classified. For example, the two sides can be jobs and machines; each machine is assigned several jobs 
but perhaps cannot take two jobs with heavy memory requirements. 

To make the problem precise, we introduce necessary notation and terminology. A set A of ap- 
plicants and a set 2 of institutes are given. Each applicant /institute has a strictly-ordered (but not 
necessarily complete) preference list over the other side. The notation indicates either strictly 
better or equal in terms of preference of an entity e € A U I while means strictly better. For 
example, if applicant a £ A strictly prefers institute i € X to another institute i' G X, we write 
i i'- The preference list of institute i is denoted as L*. The set of applicants on L* who rank higher 
(respectively lower) than some particular applicant a are written as L^^ (respectively Ll^^). 

An institute i has a capacity Q{i) G Z"*", the maximum number of applicants it can hire. It defines 
its own classification C(i) = {C^j']^'^^^\ which is a family of sets over the applicants in its preference 
list. Each class Cj G C(i) has an upperbound (?"'"(C*) G Z+ and a lowerbound g~(Cj) G Z+ U {0}, 
on the number of applicants it would hire in that class. Given a matching ^u, /x(a) is the institute 
applicant a is assigned to. We write ^{i) = (o^i, 0^2, • • • , flifc), k < Q{i) to denote the set of applicants 
institute i gets in /x, where aij are listed in decreasing order based on its preference list. In this paper, 
we will slightly abuse notation, treating an (ordered) tuple such as ^{i) as a set. 

Definition 1. Given a tuple t = (aii,ai2, ■ ■ ■ ,aife) where aij are ordered based on their decreasing 
rankings on institute i 's preference list, it is said to be a feasible tuple of institute i, or just feasible 
for short, if the following conditions hold: 

- k< Q(i); 

- given any class C] G C(i), g-(q) < |t n Cj] < g+(q). 

Definition 2. A matching ^ is feasible if all the tuples fJ.{i), i £ I are feasible. A feasible matching 
is stable if and only if there is no blocking group. A blocking group is defined as follows. Let fi(i) = 
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{cLii, ai2, - - - ,aik),k < Q{i). A feasible tuple g = Q^i2' ' ' ' ''^k-')>^ ^ k' < Q{i), forms a blocking 
group {i;g) with institute i if 

- for 1 < j < k,i ^q/^ ^(a-j) and a'-j atj; 

— either there exists 1,1 < I < k such that a'^i >-i an and i >-a'.i f^ia'u), or that k' > k. 

Informally speaking, the definition requires that for a blocking group to be formed, all involved 
applicants have to be willing to switch to, or stay with, institute i. The collection of applicants in the 
blocking group should still respect the upper and lower bounds in each class; moreover, the institute 
gets a strictly better deal (in the Pareto-optimal sense). Note that when there is no class lower bound, 
then the stable matching as defined in Definition [2] can be equivalently defined as a feasible matching 
without the conventional blocking pairs (see Lemma [T7] in Section d]) . When the class lower bound 
is present, the definition of the blocking groups captures our intuition that an institute should not 
indiscriminately replace a lower ranking applicant assigned to it with a higher applicant (with whom 
it forms a blocking pair), otherwise, the outcome for it may not be a feasible one. 

In our proofs, we often use the notation fi{i)\'^a' to denote a tuple formed by replacing a E ^(i) 
with a'. The order of the tuple )u(i)|°a' is still based on institute i's preference list. If we write n{i)\a, 
then this new tuple is obtained by adding a into fi{i) and re-ordered. In a matching fi, if a class Cj is 
fully-booked, i.e. |//(i) n Cj| = q^{Cj), we often refer to such a class as a "bottleneck" class. We also 
define an "absorption" operation: given a set B of classes, ^{B) returns the set of classes which are 
not entirely contained in other classes in B. 

Our Results It would be of interest to know how complicated the classifications of the institutes 
can be while still allowing the problem a polynomial time algorithm. In this work, we study the 
CLASSIFIED STABLE MATCHING problems whose classifications belong to a fixed set of "order types." 
The order type of a classification is the inclusion poset of all non-empty intersections of classes. We 
introduce necessary definitions to make our statement precise. 

Definition 3. The class inclusion poset P{i) = (C(i), ^) of an institute i is composed of sets of the 
elements from V: C{i) = {C\C = C] n C^, where C],Ci G C(i)]0. In P(i), C'j y cl ifC] D cl; 
and C^ ||C^ if'C^j ^ C\ and T^C^j- 

Definition 4. Let P = {Pi, P2, ■ ■ ■ ,Pk} be a set of posets. A classified stable matching instance 
{A, I) belongs to the group of P-CLASSIFIED STABLE MATCHING problems if for each poset Pj € P, 
there exists an institute i ^ Z whose class inclusion poset P{i) is isomorphic to Pj and conversely, 
every class inclusion poset P{i) is isomorphic to a poset in P. 

We call a poset a downward forest if given any element, no two of its successors are incomparable. 
Our first main result is the following dichotomy theorem. 

Theorem 5. Let P = {Pi, P2, - • • ,Pk} be a set of posets. ¥ -classified stable matching problems can be 
solved in polynomial time if every poset Pj G P is a downward forest; on the other hand, if P contains 
a poset Pj which is not a downward forest, the existence of a stable matching is NP-complete. 

^ Note that this definition allows a class to intersect itself, i.e., C — Cj nCj. This implies that C(i) 3 C(i). 
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We remark that if P is entirely composed of downward forests, then every classification C{i) must 
be a laminar famil}ll. In this case, we call the problem laminar classified stable matching 
(henceforth LCSM). 

We present an 0{m?) time algorithm for LCSM, where m is the total size of all preferences. Our 
algorithm is extended from the Gale-Shapley algorithm. Though intuitive, its correctness is difficult 
to argue due to various constraint^. Furthermore, we show that several well-known structural results 
in the hospitals/residents problem can be further generalized in LCSM. On the other hand, if 
some institute i has a classification C(z) violating laminarity, then P must contain a poset which has 
a "V" (where the "bottom" is induced by two intersecting classes in C{i) which are its parents "on 
top.") We will make use of this fact to design a gadget for our NP-complete reduction. In particular, 
in our reduction, all institutes only use upperbound constraints. Sections 2 and 3 will be devoted to 
these results. 

Our dichotomy theorem implies a certain limit on the freedom of the classifications defined by 
the institutes. For example, an institute may want to classify the applicants based on two different 
criteria simultaneously (say by research fields and gender); however, our result implies this may cause 
the problem to become intractable. 

In the second part, we study LCSM using a mathematical programming approach. Assume that 
there is no lower bound on the classes. We extend the set of linear inequalities used by Baiou and 
Balinski [3] to describe stable matchings and generalize a bin-packing algorithm of Sethuraman, Teo, 
and Qian [22] to prove that the polytope is integral. The integrality of our polytope allows us to use 
suitable objective functions to obtain various optimal stable matchings using Ellipsoid algorithm. As 
our LP has an exponential number of constraints, we also design a separation oracle. 

By studying the geometric structure of fractional stable matchings, we are able to generalize a 
theorem of Teo and Sethuraman [23]: in (one-to-one) stable marriage, given any number of stable 
matchings, if we assign every man his median choice among all women with whom he is matched in 
the given set of matchings and we do similarly for women, the outcome is still a stable matching. This 
theorem has been generalized in the context of hospitals/residents problem |5|13|22j . We prove 
that in LCSM, this theorem still holds: if we apply this "median choice operation" on all applicants, 
the outcome is still a stable matching. 

A final ramification of our polyhedral result is an answer to an open question posed by Sethuraman, 
Teo and Qian [22]: how do we describe the stable matching polytope in the classical "unclassified" 
many-to-many stable matching problem? We show this problem can be reduced to LCSM by suitable 
cloning and classifications. 

All the polyhedral results will be presented in Section 4. In Section 5 we conclude. Omitted proofs 
and details can be found in the appendix. 

1.1 Related Work 

Stable matching problems have drawn the intensive attention of researchers in various disciplines in 
the past decades since the seminal paper of Gale and Shapley For a summary, see |11I14I18] . 
Vande Vate [21] initiated the study of stable matching using mathematical programming approach; 
further developments using this approach can be found in |1I3|17|19|21|22|23] . 

^ A laminar family has no pair of intersecting classes, that is, if A,B £ JT, then either A n B = 0, or A C B, or 
B <ZA. 

^ If there is no lower bound on the classes, the proof can be significantly simplified. 

* However, this "median choice operation" on the applicants' side does not mean that all institutes will get their 
(lexicographically) median outcome. Thus, the "fairness" of the resultant matching is to a certain degree lost in our 
generalizations. See Section 4 for details. 
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Fleiner [6] studied the many-many stable matching in a much more general context. Using a fixed- 
point approach, he proved that stable matchings always exist provided that the preference of each 
entity is a substitutable choice function. Roughly speaking, such a function can be realized by imposing 
a matroid over a linear order of elements. In LCSM, supposing that there is no lower bound on the 
classes, then each laminar family is equivalent to a partition matroid. We prove that stable matchings 
always exist in this situation. Hence, our algorithm in Section 2 can be seen as a constructive proof 
of a special case of Fleiner's existence theorem. 

Abraham, Irving and Manlove introduced the student-project allocation problem [2]. It 
can be shown that in LCSM, if all classifications are just partitions over the applicants and there is 
no lower bound, our problem is equivalent to a special case of their problem. They posed the open 
question whether there is a polynomial time algorithm for their problem if there is lower bound on 
the projects (classes). Our result in Section 2 gives a partial positive answer. 

Two recent works |4|12j also consider the Hospitals/Residents problem in the context of having 
lower bounds on the hospitals' side. In [^, an interesting variation of the Hospitals/Residents problem 
which has a similar flavor to the current work is defined as follows: each hospital has its individual 
quota and sets of hospitals may also have collective quotas. 

2 An Algorithm for Laminar Classified Stable Matching 

In this section, we present a polynomial time algorithm to find a stable matching if it exists in the 
given LCSM instance, otherwise, to report that none exists. 

We pre-process our instance as follows. If applicant a is on institute i's preference list, we add a 
class C^i = {a} into C{i). Furthermore, we also add a class C| into C{i) including all applicants in 
L*. After this pre-processing, the set of classes in C{i) form a tree whose root is the CJ; moreover, 
an applicant a belongs to a sequence of classes a(C(i)) = (C*]^, C*2, ■ ■ ■ , C\^{= C|)), which forms a 
path from the leaf to the root in the tree (i.e., C*^- is a super class of C^j'' provided / < j.) For each 
non-leaf class Cj, let c(Cj) denote the set of its child classes in the tree. We can assume without loss 
of generality that g~(Cj) > Ylc- &c{c^.) 1~(.^k) non-leaf class Cj. Finally, let g'+(C|) := Q{i), 

Q'iC^) ■= Eq:ec{cj) 9~(Cfc); fo'^ applicants a G L*, g+(C*i) := 1 and i?"(C*i) := 0. 

Our algorithm finds an applicant-optimal-institute-pessimal stable matching. The applicant-optimality 
means that all applicants get the best outcome among all stable matchings; on the other hand, 
institute-pessimality means that all institutes get an outcome which is "lexicographically" the worst 
for them. To be precise, suppose that fi{i) = {aiuai2, • ■ ■ , Oifc) and fj.'{i) = (oji, aj2, • • • , dik) are the 
outcomes of two stable matchings for institute o If there exists k' < k so that Ojj = a[j, for all 
1 < J ^ fc' — 1 and aj^./ a^^,, then institute i is lexicographically better off in /_f than in fi' . 

We now sketch the high-level idea of our algorithm. We let applicants "propose" to the institutes 
from the top of their preference lists. Institutes make the decision of acceptance/rejection of the 
proposals based on certain rules (to be explained shortly). Applicants, if rejected, propose to the next 
highest-ranking institutes on their lists. The algorithm terminates when all applicants either end up 
with some institutes, or run out of their lists. Then we check whether the final outcome meets the 
upper and lower bounds of all classes. If yes, the outcome is a stable matching; if no, there is no stable 
matching in the given instance. 

How the institutes make the acceptance/rejection decisions is the core of our algorithm. Intuitively, 
when an institute gets a proposal, it should consider two things: (i) will adding this new applicant 

In LCSM, an institute always gets the same number of applicants in all stable matchings. See Theorem 1151 below. 
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violate the upper bound of some class? (ii) will adding this applicant deprive other classes of their 
necessary minimum requirement? If the answer to any of the two questions is positive, the institute 
should not just take the new applicant unconditionally; instead, it has to reject someone it currently 
has (not necessarily the new comer). 

Below we will design two invariants for all classes of an institute. Suppose that institute i gets a 
proposal from applicant a, who belongs to a sequence of classes a{C{i)) = (C*]^, C*2, • • • , CJ). We check 
this sequence of classes from the leave to the root. If adding applicant a into class C^j does not violate 
these invariants, we climb up and see if adding applicant a into C*^-^^^ violates the invariant. If we 
can reach all the way to C| without violating the invariants in any class in a(C(i)), applicant a is just 
added into institute i's new collection. If, on the other hand, adding applicant a into C*^-^-^^ violates 
the invariants, institute i rejects some applicant in C*^^_j_-^^ who is from a sequence of subclasses of 
^^^c/i can afford to lose one applicant. 

We define a deficiency number Z^(Cj) for each class S C(i). Intuitively, the deficiency number 
indicates how many more applicants are necessary for class Cj to meet the lower bound of all its 
subclasses. This intuition translates into the following invariant: 

Invariant A: A(Cj) > Ec^ec(q) ^(QO-Vq ^ C{i),c{Cj) / 0,Vi G I. 

In the beginning, A{Cj) is set to q~{Cj) and we will explain how /\(Cj) is updated shortly. Its 
main purpose is to make sure that after adding some applicants into Cj, there is still enough "space" 
for other applicants to be added into Cj so that we can satisfy the lower bound of all subclasses of 
Cj. In particular, we maintain 

Invariant B: g-(Cj) < n q| + A{C'j) < g+(Cj), VCj e C(i), Vi G I. 

We now explain how A{Cj) is updated. Under normal circumstances, we decrease A{Cj) by 1 
once we add a new applicant into Cj. However, if Invariant A is already "tight", i.e., A{Cj) = 
Sc^Gc(c*) then we add the new applicant C* without decreasing A{Cj). The same situation 

may repeat until the point that H Cj| + A{Cj) = q'^iCj) and adding another new applicant in 
Cj is about to violate Invariant B. In this case, something has to be done to ensure that Invariant B 
holds: some applicant in Cj has to be rejected, and the question is whom? 

Let us call a class a surplus class if n Cj| + A{Cj) > q~{Cj) and we define an affluent set for 
each class Cj as follows: 

$(Cj,/i(i)) = {a\a G n Cj; for each Cj, G a(C(i)) and Cj, C Cj, |^(Cj,)| + Z\(Cj,) > g"(Cj,)}. 

In words, the affluent set $(Cj,^(z)) is composed of the set of applicants currently assigned to 
institute i, part of Cj, and each of whom belonging to a sequence of surplus subclasses of Cj. In our 
algorithm, to prevent Invariant B from being violated in a non-leaf class Cj, institute i rejects the 
lowest ranking applicant a in the affluent set $(Cj,/i(i)). 

The pseudo-code of the algorithm is presented in Figure [H 
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Initialization 

0: Vi G X, vq e C{i), Z\(Cj) — g-(cj); 

Algorithm 

1: While there exists an applicant a unassigned and he has not been rejected by all institutes on his list 

2: Applicant a proposes to the highest ranking institute i to whom he has not proposed so far; 

3: Assume that a{C{i)) = (CL, C^a, ■ • • , CU= Cj)); 

4: := U {a} // Institute i accepts applicant a provisionally; 

5: For t = 2 To z / / applicant a can be added into C^i directly; 

6: If Zi(C:i) > Ec-ec(c-,) ^(^fe) Then A{C:,) ~ Aid,) - 1; 

7: If #{C',t) + A{Cl,) >q+{Cl,) Then 

8 Let %[CluK^) = {a\a G n Q^; for each C], G a(C(i)) and C], C Q^, 1^(^01 + ^(^0 > '/"(qO); 

9 Let the lowest ranking applicant in $(Cnt,/i(j)) be a^; 

10 := /i(i)\{a^} // Institute i rejects applicant a^; 
11: GOTO 1; 

12: If there exists an institute i with /i(CJ) > Then Report "There is no stable matching"; 

13: Else Return the outcome ^, which is a stable matching; 

Fig. 1. The pseudo code of the algorithm. It outputs the appUcant-optimal-institute-pessimal match- 
ing /i if it exists; otherwise, it reports that there is no stable matching. 

2.1 Correctness of the Algorithm 

In our discussion, C*^ is a class in a(C(i)), where t is the index based on the size of the class in 
a(C(i)). Assume that during the execution of the algorithm, applicant a proposes to institute i and 
when the index t of the For loop of Line 5 becomes I and results in being rejected, we say applicant 
a is stopped at class C*^, and class C*^ causes applicant to be rejected. 
The first lemma describes some basic behavior of our algorithm. 

Lemma 6. (i) Immediately before the end of the while loop, Invariants A and B hold, 
(a) Let applicant a be the new proposer and assume he is stopped at class C*^. Then 

(iia) Between the time interval that he makes the new proposal and he is stopped at C*;, Z\(C*j) 
remains unchanged, for all 1 < t < I; moreover, given any class C^f, 2 < t < I, Z\(C*^) = 

(iib) When a is stopped at a non-leaf class C'^^i, ${C^i, fi{i)) ^ 0; in particular, any class C*j, 
1 < t < I — 1, is temporarily a surplus class. 

(iii) Immediately before the end of the while loop, if class Cj is a non-leaf surplus class, then ^(Cj) = 
Sc^Gc(Cj) ^(Cfc)- 

(iv) Suppose that applicant a is the new proposer and C^^i G a(C(i)) causes applicant to be rejected 
and a^{C{i)) = {C'^^fi, Cat2' " " " ' ^atit(~ ^ai)^ ' ' ')■ Then immediately before the end of the while 
loop, A{Cl,^,) = f:ciec(C' )^iCl), for all 2 < t' < l^ ; moreover, \^i{i) n C^^^tl + ^(C^lt^t) = 

Proof, (i) can be proved by induction on the number of proposals institute i gets. For (iia), since 
Invariant A is maintained, if Z\(C*j) is decreased for some class C*^, 1 < t < I, the algorithm will 
ensure that applicant a would not be stopped in any class, leading to a contradiction. Now by (iia), 
the set of classes {C*^}'~| are (temporarily) surplus classes when applicant o is stopped at C*;, so 
$(C*;, /x(i)) 7^ 0, establishing (iib). Note that this also guarantees that the proposed algorithm is never 
"stuck." 
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(iii) can be proved inductively on the number of proposals that institute i gets. Assuming a is the 
new proposer, there are two cases: (1) Suppose that applicant a is not stopped in any class. Then a 
class C*^ G a(C(i)) can become surplus only if the stated condition holds ; (2) Suppose that applicant 
a is stopped in some class, which causes to be rejected. Let the smallest class containing both a and 
be C*;,. Applying (iia) and observing the algorithm, it can be verified that only a class C C*^ 
can become a surplus class and for such a class, the stated condition holds. 

Finally, for the first part of (iv), let C*;, denote the smallest class containing both a and . Given 
a class C*tt" — — ^a/' gives the proof. If C*|^, C C*;,, observe that the former must 

have been a surplus class right before applicant a made the new proposal. Moreover, before applicant 
a proposed, (iii) implies that for a non-leaf class C*|^, C C*;,, the stated condition regarding the 
deficiency numbers is true. The last statement of (iv) is by the algorithm and Invariant B. □ 

Lemma 7. Assume that a\'C{i)) = (C^ti' ^at2' " " " '^atit'''')- During the execution of the algo- 
rithm, suppose that class C*|^| causes applicant to he rejected. In the subsequent execution of the 
algorithm, assuming that iJ,{i) is the assignment of institute i at the end of the while loop, then 
there exists l^, where l^ > l^ such that \fJ.{i) H C^^l + ^(^at;t) ~ ^^IC^tit)'' furthermore, for 
all 2 < t < l^, all applicants in $(C*|j, /x(i)) rank higher than a"^ . Moreover, for all 2 < t < l^. 

Proof. We prove based on the induction on the number of proposals institute i receives after is 
rejected. The base case is when is just rejected. Let = l^ . Then it is obvious that all applicants 
in the affluent sets $(C*|^, /i(i)), 2 < t < l^, rank higher than and the rest of the lemma holds by 
Lemma [6|^iv) . 

For the induction step, let a be the new proposer. There are four cases. Except the second case, 
we let l^ remain unchanged after a's proposal. 

— Suppose that a C'^t^t ^^^^ does not cause anyone in C*|^j to be rejected. Then the proof is 
trivial. 

— Suppose that a C'*tit stopped in class C*;, which causes an applicant a* € C^^-tn to 
be rejected, a* must be part of the affluent set $(C*|^(, ;u(z)) before a proposed. By induction 
hypothesis, a* >-i a^. Moreover, since a* is chosen to be rejected, all the applicants in the (new) 
affluent sets $((7*^^, /x(z)), for each class C*|^, where C^^^j C (7*^^ C C*;, rank higher than a*, 
hence, also higher than a^. Now let C*; be the new C^^x ^^'^ ^^^^ °^ lemma follows from 
Lemma O^iv) . 

Suppose that a G C^^Ut stopped in C^^^j or any of its subclasses. We argue that a must 

be accepted without causing anyone to be rejected; moreover, the applicants in all affluent sets 
$(C*|^, //(«)), for all 1 < t < /-I- remain unchanged. Let the smallest class in a^{C{i)) containing a be 
^atf before a proposed, the induction hypothesis states that |/i(i) nC*|^j| + Z\(C*^^j) = 

g"'"(C*|^(). As applicant a is not stopped at C^^^j, the set of values Z\(C*|J, I < t < l^, must have 
decreased during his proposal and this implies that he will not be stopped in any class. 
Now let a(C(.)) = (Q^, • ■ ■ , C^,, q(,+,)(= C^,^-), • • • )• Since Z\(q,^-) = Ec^ec(c^^-) before 

applicant a proposed by the induction hypothesis, for Z\(C*^-) to decrease, Z^(C*;) must have 
decreased as well. Choose the smallest class C*;« C C*^- whose value Z\(C*;,) has decreased during 
a's proposal. We claim that C*;, must have been a non-surplus class before and after applicant 
a's proposal. If the claim is true, then all the affluent sets ${C^^^, fi{i)), for all 1 < t < l^, remain 
unchanged after applicant a's proposal. 
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It is obvious that C*^, 7^ C*]^. So assume that C*^, is a non-leaf class. Suppose for a contra- 
diction that C*;* was a surplus class before a proposed. Lemma [D^iii) implies that Z\(C*|^, ) = 
X]c'ec(C' ) ^(^fc) before a proposed. Then for Z\(C*|j,) to decrease during a's proposal, Z\(C*|^^,_-|^^^ 

must have decreased as well. But then this contradicts our choice of C%,.. So we establish that 
^ai* ^'-'^ surplus and remains so after a's proposal. 
— Suppose that a € C^^j^f and when he reaches a subclass of C*|^j or the class itself, the latter causes 
some applicant a* to be rejected. To avoid trivialities, assume a ^ a*. Let the smallest class in 
a^(C(i)) containing a be C^^^- and the smallest class in a^(C(z)) containing a* be C^-\i, - Below we 
only argue that the case that C*^- C C*|^,. The other case that C*|j, C C*^- follows essentially 
the same argument. 

After a's proposal, observe that only the affluent sets $(C*|^, I < t < I*, can have new 
members (who are from the child class of C*^- containing a). Without loss of generality, let G be 
the set of new members added into one of the any above sets. 

To complete the proof, we need to show that either G = or all members in G rank higher than 
a^. If before applicant a proposed, a* belonged to a sequence of surplus classes C*,^ C C*ti*' 
was also part of the affluent set $(C*|j, , /i(i)) or part of nC*|-^ before a proposed. By induction 
hypothesis, a* a) . Observing LemmalD^iib), all applicants in G must rank higher than a*, hence 
also than . On the other hand, if a* belongs to some class C*.^ C C^^^^ which was not surplus 
before a proposed, then C*^- = C*.;* and C**^ must also contain a and remain a non-surplus class 
after a's proposal. In this case G = %. □ 

The following lemma is an abstraction of several counting arguments that we will use afterwards. 

Lemma 8. Let each class Cj be associated with two numbers a^- and and q~{Cj) < Oi^j,Pj < 
q'^{Cj). Given any non-leaf class C*, a* = J2ciec{C^.) ^^'^ (^j — X]c^ec{C*) i^L' moreover, if Pj = 
Sc'6c(C') (^k then such a non-leaf class Gj is said to be tight in j3. If (3j > q~{Cj), then Cj has to be 
tight in j3. 

(i) Given a non-leaf class C^^^^ with a^t^t < l^liHV '"^^ can find a sequence of classes C^^^i D • • • D C*ti' 
where a%, < /?%,, for 1 < t < /t. 

ait ' aif ■' — — . . . . 

(a) Given a non-leaf class C* with a]. < /?* , suppose that there exists a leaf class C^^-j^ C C* such that 
^ ^\<i'Y Moreover, all classes C*^^ are tight in j3, where C*^-^ C C*^^ C C^, then we can find 
a class C*/, where G^^^-^ C C*/ C C*, a]., < (3^.,, and two sequences of classes with the following 
properties: 

(ua) Gl,^ CGI,^C---C Gl,^, C Gl„ where a\,^ > forl<t< if; 
(iib) Gi, D Gl,i, D---DGi,^, where a\,^ < [3^,^, for I < t < ll 



Proof For (i), since q (G^^^^) < a^^^^^ < class C'^^t is tight in (3. Therefore, Ec^gc(C'^^^) «fc = 

"at/t < /^at/t = Hci^cici^^^) Pk- By counting, there exists a class C*t(;t-i) ^ c(C*t;t) with Q"(C^t(;t_i)) < 
'^at(it-i) ^ '^at(/t-i)- Repeating the same argument gives us the sequence of classes. 

For (ii), let us climb up the tree from G^^^^ until we meet a class C*, C with q^, < This 
gives us the sequence of classes stated in (iia). 

Now since the class C*, is tight in /?, Y^ci&c{c\)^k ~ "^x' — (^x' ~ J2c'Igc{c\) f^l- Moreover, as 
C*^;^ € c(C*,) and a^^^^ > by counting, we can find another class C^^^t G c(C*,)\{C*^^^} such 

that > a^t^t — 9~(^at/t)- Now applying (i) gives us the sequence of classes in (iib). □ 
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We say that {i; a) is a stable pair if there exists any stable matching in which appHcant is assigned 
to institute i. A stable pair is by-passed if institute i rejects applicant a during the execution of our 
algorithm. 

Lemma 9. During the execution of the algorithm, if an applicant a"^ is rejected by institute i, then 
{i;a'^) is not a stable pair. 

Proof. We prove by contradiction. Assume that {i;a'^) is the first by-passed stable pair and there 
exists a stable matching in which ij,'f'{a'^) = i. For each class Cj G C(i), we associate two numbers 
a* := \fi'^{i) n Cjl and /?j := \fi{i) Pi Cj| + Z^(Cj). Here A{-)s are the values recorded in the algorithm 
right after a'^ is rejected (before the end of the while loop); similarly, fi{i) is the assignment of i at 
that point. 

It is obvious that a^^-^ > P^^^^ and the class C* causing a'^ to be rejected is not C*^-^. By 
Lemma El^iv) , all classes C*^^ are tight in /?, where C^^j^ C C*,^^ C C*. It can be checked all the 
conditions as stated in Lemma [8)^ii) are satisfied. In particular, /?* = q~^{Cl) > moreover, if 

> '?~(C'])> C*] must be tight (by Lemma [H in)). 

So, we can find two sequences of classes {C^^jlJ^^;^ and {C*|j}'L|, where C**^^,^, C*|^| G c(C*/) and 
C*, C C*, with the following properties: 

g+(q,,) > l/^nO n > n + zl(q,J > g-(q,J, Vi, \<t<l<i>- 
q-{ci,,) < |/(i) n c^tj < |//(i) n c^ttl + ^Ku) < 9+(c:tt), vt, i < i < /t. 

The second set of inequalities implies that the classes {C*^^}jL| are surplus in /i. Thus there exists 
an applicant € {iJ.{i)\fi'^ (i)) nC*^-^. Since {i]a'^) is the first by-passed stable pair, i >-^\ fi'^{a^) and 
since a"^ is rejected instead of a^, >-i a'^. Now observe the tuple fi'^ (i)]'^''' is feasible due to the 
above two sets of strict inequalities. Thus we have a group (z;/i<^(i)|'^ at) to block fi^, a contradiction. 

□ 

Lemma 10. At the termination of the algorithm, if there exists an institute i such that Z^(CJ) > 
0, there is no stable matching in the given instance. 

Proof. Suppose, for a contradiction, that there exists an institute i with A{C'^) > and there is a 
stable matching n^. Let be the assignment when the algorithm terminates. By Lemma [9l if an 
applicant is unmatched in ^, he cannot be assigned in fi'^ either. So < In the following, /i(-)s 
refer to values recorded in the final outcome of the algorithm. Consider two cases. 

— Suppose that |/u'^(«)| > |/-f(i) H CJ|. Then as l/i*^! < we can find another institute i' ^ i 
such that |/u'^(i')| < |/u(i') n q'|. For each class Cf £ C(z'), let a| := |/i*(z') n Cf\ and := 
\fi{i')ncj\ + A{C'j). It can be checked that the condition stated in Lemma El^i) is satisfied (note 
that those fulfill the condition due to Lemma [6l^iii)). Therefore, we can find a sequence of 
classes {C^^^}^^i, where C^'^^i = C|', and 

|/(.') n cftj < n cftj + z\(cft,) < g+(c:;,), vt, i<t<i\ 

where the second inequality follows from Invariant B. Then there exists an applicant G 
(/i(i')V'^(«')) n Cfti- By Lemma El i' ^^t M'^(at), giving us a group (i'; /i<^(i')|a^) to block n'i', 
a contradiction. Note the feasibility of fj,'^{i')\a^ is due to the above set of strict inequalities. 
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— Suppose that < \fJ-{i) n CJ|. We first claim that CJ must be a surplus class in fi{i). If not, 

then g-(q) = z\(q) + \fi{i) n cjl > n q|, implying that > g-(q) > n qi, a 

contradiction. So q is a surplus class, and by Lemma [U^iii) , 

E l/«nc^|<|Mi)nq|<IMOnqi + z^(q)= J] IM^) nc^i + 

For each class q G C(i), let a* := |^<^(i)nq| and /3j := |^(i)nq|+Z^(q) and invoke LemmaEJi). 
The above inequality implies that aj < /3J and note that by Lemma [6l^iii) , the condition regarding 
(3 is satisfied. Thus we have a sequence of surplus classes qi^|(= q) D ■ ■ ■ D q^-^ so that 

^-(qti) < n < 1^(0 n qtti + ^Ku) < Q^ic:u)yt, i<t<i\ 

implying that there exists an applicant € n q^-^ and i :^^t fJ-'^{a^) by virtue 

of Lemma [H The tuple fi'^{i)\a^ is feasible because of the above set of strict inequalities. Now 
{i; fj.'^ {i)\a'^) blocks fi'^, a contradiction. □ 

Lemma 11. Suppose that in the final outcome fi, for each institute i ^I, ^(q) = 0. Then fi is a 
stable matching. 

Proof. For a contradiction, assume that a group {i]g) blocks /u. Let a'^ to be the highest ranking 
applicant in g\fi{i). Since a"^ is part of the blocking group, he must have proposed to and been 
rejected by institute i during the execution of the algorithm, thus i >-^g> fJ'io-'^)- By Lemma [71 there 
exists a class q^^j such that n q^^J + -^(q^^t) = n q^^J = g+(q^^(). Moreover, it is 
obvious that \g fl q,^]^| > n C^^^^\. We now make use of Lemma O^ii) by letting a* := |(7 n q| 
and f3j := \n{i) n q| for each class q E C(i). Note that all classes are tight in /?, q^^^ C q,^,;}> and 
|/x(i)nq^^f I = (?"'"(q^^() > l^nq^^j |, satisfying all the necessary conditions. Thus, we can discover a 

sequence of classes {q^^}'^]^ stated in Lemma [8)^iib), where qt;t ^ c{C^^^^) and q,/,-^ C q^^ C q^^j, 
such that 

g-(qti) < b n qtti < 1^(0 n qtti < 9+(C^:tt), Vj, 1 < t < /t, 

and there exists an applicant a) G {fi{i)\g) fl qtj^- The above set of strict inequalities mean that all 
classes qtjj ^ ^ t < , are surplus classes in fi. Then forms part of the affluent set $(q^^, 
By Lemma [71 they all rank higher than a*^. This contradicts our assumption that a"^ is the highest- 
ranking applicant in g\n{i). □ 

Lemma 12. Suppose that in the final outcome ^, for each institute i ^I, ^(q) = 0. Then ^ is an 
institute-pessimal stable matching. 

Proof. Suppose, for a contradiction, that there exists a stable matching ^jl'^' such that there exists 
an institute i which is lexicographically better off in ^ than in . Let be the highest ranking 
applicant in ^{i)\^i't> {i) . By Lemma [3 i ;-^t /""^(i)- If n qt J < \^i{i) n qt J < q+{Ci^^), for 

all classes q^^ S a"l'(C(i)), then (i; /Lt'^(i)|a'^) blocks /u"^, a contradiction. So choose the smallest class 
q G aUC(«)l such that n q| > n q|. It is clear that q D q^^. 

Now we apply Lemma [Hj^ii) by letting q*- := \^{i) fl q| and /?* := H q| for each class 

q G C(z). It can be checked all conditions stated in Lemma [S^ii) are satisfied. So there exists a class 
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C*, such that C*|-^ C C*, C and we can find two sequences of classes {C^^^^}[^i and {C'*tt}' 
where 0"^^^^^^ ^atn ^ c(C*./), with the following properties: 




The second set of inequalities implies that we can find an applicant a'^ G n C^^-^- 

Recall that we choose to be the highest ranking applicant in (i) , so >-{ a'^. Now we have 

a group (i; iJ,'^{i)\°'''' a'^) to block fi'f' to get a contradiction. The feasibility of fi'^ (i)]'^''' is due to the 



Based on Lemmas [9l [101 fill ^^id [121 we can draw the conclusion in this section. 

Theorem 13. In 0{m?) time, where m is the total size of all preferences, the proposed algorithm 
discovers the applicant-optimal-institute-pessimal stable matching if stable matchings exist in the given 
LCSM instance; otherwise, it correctly reports that there is no stable matching. Moreover, if there is 
no lower bound on the classes, there always exists a stable matching. 

To see the complexity, first note that there can be only 0{m) proposals. The critical thing in 
the implementation of our algorithm is to find out the lowest ranking applicant in each affluent set 
efficiently. This can be done by remembering the lowest ranking applicant in each class and this 
information can be updated in each proposal in 0{m) time, since the number of classes of each 
institute is 0{m), given that the classes form a laminar family. 

2.2 Structures of Laminar Classified Stable Matching 

Recall that we define the "absorption" operation as follows. Given a family of classes B, ^(B) returns 
the set of classes which are not entirely contained in other classes in B. Note that in LCSM, 3?(-B) 
will be composed of a pairwise disjoint set of classes. 

We review the well-known rural hospitals theorem |8|15j . 

Theorem 14. (Rural Hospitals Theorem) In the hospitals/residents problem, the following holds. 

(i) A hospital gets the same number of residents in all stable matchings, and as a result, all stable 
matchings are of the same cardinality. 

(ii) A resident who is assigned in one stable matching gets assigned in all other stable matchings; 
conversely, an unassigned resident in a stable matching remains unassigned in all other stable 
matchings. 

(Hi) An under- subscribed hospital gets the same set of residents in all other stable matchings. 

It turns out that rural hospitals theorem can be generalized in LCSM. On the other hand, if 
some institutes use intersecting classes in their classifications, rural hospitals theorem fails (stable 
matching size may differ). See the appendix for such an example. 

Theorem 15. (Generalized Rural Hospitals Theorem in LCSM) Let ^ be a stable matching. Given 
any institute i, suppose that B is the set of bottleneck classes in and D is the subset of classes 
in C{i) such that ^(B) U D partitions L*. The following holds. 



above two sets of strict inequalities. 



□ 
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(i) An institute gets the same number of applicants in all stable matchings, and as a result, all 

stable matchings are of the same cardinality, 
(a) An applicant who is assigned in one stable matching gets assigned in all other stable matchings; 

conversely, an unassigned applicant in a stable matching remains unassigned in all other stable 

matchings. 

(Hi) Every class CI. € ^{B) U D has the same number of applicants in all stable matchings. 

(iv) In a class CI, C C (z D, or in a class CI. which contains only classes in D, the same set of 

applicant in class C^ will be assigned to institute i in all stable matchings. 
(v) A class CI can have different sets of applicants in different stable matchings only if C\ ^ C ^ 

^{B) orCl DCs ^{B). 



Proof. We choose /x^ to be the apphcant-optimal stable matching. 

Claim A: Suppose that a € fi^i)\fi{i). Then there exists a class C*; G a(C(i)) such that (i) 
W) nCli\= q+{Cl,), and (ii) a G C^, C C G ^{B). 

Proof of Claim A. If for all classes C*^ G a(C(i)), nC*j| < g^(C*^), then as ^'^ is applicant- 
optimal, i l^{o)^ so («;/u(i)|a) blocks /i, a contradiction. This establishes (i).(ii) follows easily. □ 

Let B (1 B he the subset of these bottleneck classes containing at least one applicant ij){i)\^{i). 
By Claim A(n), ^{B) C ^{B). This implies that for all classes Cl G {?R:{B)\^{B)) U D , \fi{i)nCl\ > 
n C|.|. Combining this fact with Claim ^(ii), we have 

> Y l/^^Wnc^l+ Y ^^^^al) (*) 
Y \fi\i)nci\+ Y I^^WnQI 

C^G{R{-B)\K{-B))UD C^gR{-B) 
= |^t(,)|. 

Thus, > and it cannot happen that > otherwise, there exists an applicant who 
is assigned in fi but not in fi^ . This contradicts the assumption that the latter is applicant-optimal. 
This completes the proof of (i) and (ii) of the theorem. 

Since |^| = Inequality (*) holds with equality. We make two observations here. 

Observation 1: For each class Cl. G ^{B), it is also a bottleneck in fi^{i). 
Observation 2: an applicant a G fi"^ {i)\fi{i) must belong to a bottleneck class in fJ,^{i). 

Let B^ be the set of bottleneck classes in /U^(i) and choose so that ^{B^) U partitions L*. 
By Observation 2, each applicant in //^(i) nC^, where C^ G D\ must be part of So for each class 
Cl G D\ \n{i) n Cl\ > \n^{i) n Ci\. We claim that it cannot happen that \n{i) n > Ifi^i) n Cl\. 
Suppose not. Then by (i), there are two possible cases. 

— There exists another class C|,, G so that |/Lt(i)nC^j I < |/it(i)nC|,|. Then we have a contradiction 
to Observation 2. 
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- There exists another class C^, G ^{B^) so that \fi{i)r\Cl, \ < nC|,, |. For each class Cj € C(«), 
let a] := \fi{i)nCj\ and Z?* := |/it(i) nq|. Then we can invoke Lemma [8]^i) and find an applicant 
a't' G n^{i)\iJ.{i) so that for each class C^^^ G a'^(C(i)), C^^^ C C|,, n C^^J < n C*^J < 
g"'"(C*^^). Then by Claim A(ii) and Observation 1, there must exist another class C^„ G 3ft(i?) 
containing a"*^ and C^„ D C^,. By Observation 1, C^„ is also a bottleneck class in //^(i). This 
contradicts the assumption that CI, G 

So we have that for each class CI G , n = H C^J. For each class C^, G i?^, we can 

use the same argument to show that n Cl\ = |/U^(i) n C^,|. This gives us (iii) and (iv). (v) is a 
consequence of (iv). □ 

3 NP-completeness of P-Classified Stable Matching 

Theorem 16. Suppose that the set of posets P = {Pi, P2, ■ ■ ■ ,Pk} contains a poset which is not a 
downward forest. Then it is NP-complete to decide the existence of a stable matching in V-classified 
stable matching. This NP-completeness holds even if there is no lower bound on the classes. 

Our reduction is from ONE-IN-THREE SAT. It is involved and technical, so we just highlight the 
idea here. As P must contain a poset that has a 'V in it, some institutes use intersecting classes. In 
this case, even if there is no lower bound on the classes, it is possible that the given instance disallows 
any stable matching. We make use of this fact to design a special gadget. The main technical difficulty 
of our reduction lies in that in the most strict case, we can use at most two classes in each institute's 
classification. 

4 Polyhedral Approach 

In this section, we take a polyhedral approach to studying LCSM. We make the simplifying assump- 
tion that there is no lower bound. In this scenario, we can use a simpler definition to define a stable 
matching. 

Lemma 17. In LCSM, if there is no lower bound, i.e., given any class Cj, q~{Cj) = 0, then a stable 
matching as defined in Definition\Mcan be equivalently defined as follows. A feasible matching fi is sta- 
ble if and only if there is no blocking pair. A pair {i, a) is blocking, given that = (aji, 0^2, • • • , flifc); 
k < Q{i), if 

- i >-a ^(a); 

- for any class C^ G a{C{i)), 11^1^0 fi{i) H C'J < q+iCl). 

The definition of blocking pairs suggests a generalization of the comb used by Baibu and Balin- 
ski [3]. 

Definition 18. Let P = I x A denote the set of acceptable institute- applicant pairs. The shaft S{A'^), 
based on a feasible tuple A'' of institute i, is defined as follows: 

S{A') = {{i, a')er: VCj G a'(C(i)), |L!_^, nA'nCj\< q^{C})}. 
The tooth T{i, a) is defined for every (i, a) (z P as follows: 

T{i,a) = {{i',a) eP:i'ha i}- 
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In words, {i,a') forms part of the shaft S{A^), only if the collection of a' and all applicants in 
ranking strictly higher than a' does not violate the quota of any class in a'(C(i)). We often refer to 
an applicant a G ^* as a tooth- applicant. 

We associate a |-r|-vector (or simply x when the context is clear) with a matching fi: xf^ = 1 
if fi{a) = i, otherwise, x^^ = 0. Suppose that F C F. Then x{r) = a)&r^i'^- define a comh 
K{i,S{A^)) as the union of the teeth {T{i,ai)}^^^j^i and the shaft S{A^). 

Lemma 19. Every stable matching solution x satisfies the comb inequality for any comb K{i, S{A^)): 
x{K{i,S{A')) = x{S{A')) + x{T{i,aj)\{i,aj}) > \A'\. 



It takes a somehow involved counting argument to prove this lemma. Here is the intuition about 
why the comb inequality captures the stability condition of a matching. The value of the tooth 
x{T{i,a)) reflects the "happiness" of the applicant a € A^. li x{T{i,a)) = 0, applicant a has reason 
to shift to institute i; on the other hand, the values collected from the shaft x{S{A'')) indicates the 
"happiness" of institute i: whether it is getting enough high ranking applicants (of the "right" class). 
An overall small comb value x{K{i, S{A^))) thus expresses the likelihood of a blocking group including 
i and some of the applicants in A^. 

Now let denote the set all combs of institute i. We write down the linear program: 

Xia<l,yaeA (1) 

i:{i,a)er 

Y xia<q+ic})yiei,yqecii) (2) 

a:{i,a)er,aGC^. 

x{K{i,S{A'))) = Y Xia>\A'\yK{i,S{A'))eK^yieI (3) 

{i,a)eK{i,S(A^)) 

Xia>0,y{i,a) e r (4) 

Suppose there is no classification, i.e., Hospitals/Residents problem. Then this LP reduces to 
the one formulated by Baiou and Balinski [3]. However, it turns out that this polytope is not integral. 
The example in Figure [2] demonstrates the non-integrality of the polytope. In particular, observe that 
since fi is applicant-optimal, in all other stable matchings, applicant 03 can only be matched to is. 
However, the value Xj^aa = 0.2 > indicates that x is outside of the convex hull of integral stable 
matchings. 

Here we make a critical observation. Suppose that in a certain matching /x"^, applicant 03 is 
assigned to ii. Then 02 cannot be assigned to ii due to the bound q'^{Cl) (see Constraint ([2|)). If fi'f' 
is to be stable, then 02 must be assigned to some institute ranking higher than ii on his list (in this 
example there is none), otherwise, (i, //'^(ii)|"'^a2) is bound to be a blocking group in ^'f'. Thus, the 
required constraint to avoid this particular counter-example can be written as 

x{T{ii,a2)\{ii,a2}) > Xi^aa- 

We now formalize the above observation. Given any class Cj E C(i), we define a class-tuple 
t) = (oji, «i25 ■ ■ ■ , o-iq+(c^.))- Such a tuple fulfills the following two conditions: 
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Institute Preferences Classifications Class bounds 

ii:aia6a7a2a3 Cl = {a2,a3} (5(ii) = 2, q+(Ci) = f 

i2M4a7 Q(i2) = 1 

13:0204 Q{h) = 1 

14:0506 Q(i4) = 1 

isMsasarai Ci= {03,05} Q(i5) = 2, q+(Cf) = f 



Applicant Preferences 

01:2521 

02:1113 

03:1115 

04:1312 

05:1514 

06:2421 

07:121125 



(Applicant-optimal stable matching) /i — {{ii; 0.2, a^), (12; 07), (13; 04), (24; 05), (25; oi, 03)} 
(A fractional matching x that is not inside the convex hull of integral stable matchings) 

^iiai — 0.5, Xi-^aQ — 0.8, Xi-^ay — 0.2, Xi-^^a2 — 0.3, Xi-^^a^ — 0.2. 
^22^4 — O.T, Xi2aj — 0.3. 
•^i^a2 — O.T, Xir^a^^ — 0.3. 
^i^a^ — 0.8, Xi^aQ — 0.2. 

^i^a^ — 0.8, Xic^a^ — 0.2, Xij^ay — 0.5, Xi^ai — 0.5. 

Fig. 2. An example showing that the polytope formed by Constraints is not integral. Since /i 

is applicant-optimal, in all other stable matchings, applicant 03 can only be matched to i^. However, 
the value Xj^as = 0.2 > indicates that x is outside of the convex hull of integral stable matchings. 



2. if Cj is a non-leaf class, then given any subclass CI of Cj, |t* n < q~^{Cl). 

Let L^^i denote the set of applicants ranking lower than all applicants in t*- and the set of 
applicants ranking at least as high as the lowest ranking applicant in t*-. 

Lemma 20. Every stable matching solution x satisfies the following inequality for any class-tuple t*-.' 



^ x{T{i,aij)\{i,aij}) > ^ 



As before, it takes a somehow involved counting argument to prove the lemma but its basic idea 
is already portrayed in the above example. Now let T* denote the set of class-tuples in class G C(i) 
and L* i denote the set of applicants ranking lower than all applicants in t). We add the following 
sets of constraints. 



^ x(T{i,aij)\{i,aij}) > 



E 

1 -i.t'-. 



(5) 
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Let Pfsm denote the set of all solutions satisfying ([I])-© and Pgm the convex hull of all (integral) 
stable matchings. In this section, our main result is Pfsm = Psm- We say {i,a) are matched under x 
if Xia > 0. 

Definition 21. Let x € Pfsm. md fii{x) he the set of applicants that are matched to institute i under 
X. Let fii{x) he composed of aii,ai2, ■ ■ ■ , ordered hased on the decreasing preference of institute i. 

1. Define Hi{x) as a tuple composed of applicants chosen hased on the following procedure: adding 
Oij greedily unless adding the next applicant into Hi{x) will cause Hi{x) to violate the quota 
of some class. Equivalently, an ^ Hi{x) only if there exists a class Cj G aj;(C(i)) such that 
mx)n{au}[z\\=q+{q). 

2. Define Ei{x) as a tuple composed of applicants for whom institute i is the most preferred institute 
that they are matched under x, i.e., an applicant a G Ei{x), if x{T{i, a)\{{i, a)}) = and Xia > 0. 
The order of the applicants in Ei{x) is hased on the decreasing preference of institute i. 

Lemma 22. Ei{x) is feasihle for institute i. 

Proof. We need to show that given any class C* G C(i), \Ei{x)riCj\ < q~^{Cj). We proceed by induction 
on the height of Cj in the tree structure of C(i). The base case is a leaf class. If |£'j(x)nC*| > q~^{Cj), 
form a class-tuple by picking the first q~^{Cj) applicants in £'j(x)nCj. Then Constraint ([5]) is violated 
in such a class-tuple. For the induction step, if \Ei{x)riCj\ > q'^{Cj), again choose the q~^{Cj) highest- 
ranking applicants in Ei{x) n C* and we claim they form a class-tuple of Cj, the reason being that by 
induction hypothesis, given any C Cj, \Ei{x)(~\C1\ < g~^(C^,). Now Constraint ([5]) is again violated 
in such a class-tuple. □ 

Lemma 23. Suppose that x G Pfsm- 

(i) For each institute i £2, we can find two sets U and V of pairwise disjoint classes so that U UV 
partitions L* and all applicants in Qi{x)\Hi{x) helong to the classes in U . Moreover, 

(la) \H,{x)\ = Ecieu ^^i^^l) + Y^ci^v ^x) n C^|; 

(ih) for each class CI G U, \Hi{x) H = \Ei{x) D Cl\ = q~^{Cl); for each class CI £ V and each 

applicant a G CI, if Xia > 0, then xia = 1; 
(ic) for each class CI G U, Y^aaci ^ia = q^{Cl). 

(a) For every applicant a G Hi{x), x{T{i,a)) = '^i^j^ia = l/ moreover, given any two institutes i, 

i' Hi{x)nHi>{x) = 0. 
(Hi) \Hi{x)\ = \Ei{x)\ for all institutes i 
(iv) X^dg^^a^ia = |-E'i(a;)| for all institutes i £l. 

Proof. For (i), given any applicant a G f2i{x)\Hi{x) , by Definition [211 there exists some class G 
a{C{i)) for which \Hi{x) n C*| = q~^{Cj). Let B be the set of classes Cj which contain at least one 
applicant in Qi{x)\Hi{x) and |Cj n Hi{x)\ = q~^{C!j). Let U := ?fi{B) and choose V in such a way so 
that U UV partitions L*. Now (ia) is a consequence of counting. We will prove (ib)(ic) afterwards. 

For (ii), by definition of Hi(x), none of the applicants in Qi{x)\Hi{x) contributes to the shaft 
x{S{Hi{x))). As a result, for Constraint ([3|) to hold for the comb K{i, S{Hi{x))), every tooth-applicant 
a G Hi[x) must contribute at least 1, and indeed, by Constraint ([T]), exactly 1. So we have the first 
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statement of (ii). The second statement holds because it cannot happen that x{T{i, a)) = x{T{i', a)) = 
1, given that Xia > and Xj/^ > 0. 

For (iii), By Definition 12 H all sets Ei[x) are disjoint; thus, every applicant who is matched under 
X belongs to exactly one Ei{x) and at most one Hi{x) by (ii). Therefore, Yli^i ^ Siex 

and we just need to show that for each institute i, < \Hi{x)\, and this follows by using (ia): 

1^.(^)1= E E l^.(^)ncil> E mx)nci\ + E \E,{x)nci\ = \Ei{x)l (6) 

where the inequality follows from Lemma [22] and the fact all applicants in i7i{x)\Hi{x) are in classes 
in U. So this establishes (iii). Moreover, as Inequality ([6]) must hold with equality throughout, for 
each class CI G V, if applicant a £ CI is matched to institute i under x, he must belong to both 
Hi{x) and Ei{x), implying Xia = 1; given any class CI € U, \Hi{x) nC^| = \Ei{x) nC|,| = q'^{Cl). So 
we have (ib). 

For (iv), consider the comb K{i, S{Ei{x))). By definition, x{T{i, a)\{{i, a)}) = for each applicant 
a G Ei{x). So 

x{K{i,S{E,{x))))=x{S{Ei{x))) 

= E \Ei{x)r\ci\ + E E 

Cl&V CleU a'eCl,ii,a')eS{E,{x)) 

< E mx)nci\+ E q'^{ci) = mx)\, 

where the inequality follows from Constraint ([2]) and the rest can be deduced from (ib). By Con- 
straint (l3|), the above inequality must hold with equality. So for each class € U, X^^'ecj. (i a')eS{Ei{x)) ~ 
Sa'GC^ Xja' = Q^iCl), giving us (ic) and implying that there is no applicant in € [/ who is matched 
to institute i under x ranking lower than all applicants in Ei{x) fl C^. The proof of (iv) follows by 

E^-= E E^-+ E E^-= E mx)nci\+ e q+{ci) = mx)\. 

aeA Cl&Va&Cl Cl^Ua&Ci C^eV C^GC/ 

□ 

Packing Algorithm 

We now introduce a packing algorithm to establish the integrality of the polytope. Our algorithm is 
generalized from that proposed by Sethuraman, Teo, and Qian [22]. Given x € Pfsm, for each institute 
i, we create |-E^i(3;)| "bins," each of size (height) 1; each bin is indexed by (i, j), where 1 < j < \Ei{x)\. 
Each Xia > is an "item" to be packed into the bins. Bins are filled from the bottom to the top. 
When the context is clear, we often refer to those items Xia as simply applicants; if applicant a G Cj, 
then the item Xia is said to belong to the class C*. 

In Phase 0, each institute i puts the items Xia, if a € Hi{x), into each of its |£^j(x)| bins. In the 
following phase, t = 1, 2, • • • , our algorithm proceeds by 

— first finding out the set Lj of bins with maximum available space; 
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— then assigning each of the bins in Lt one item. 

The assignment in each phase proceeds by steps, indexed by / = 1, 2, • • • ■,\Lt\. The order of the 
bins in Lt to be examined does not matter. How the institute i chooses the items to be put into its 
bins is the crucial part in which our algorithm differs from that of Sethuraman, Teo, and Qian. We 
maintain the following invariant. 

Invariant C: The collection of the least preferred items in the |£'i(rE)| bins (e.g., the items 
currently on top of institute i's bins) should respect of the quotas of the classes in C(i). 

Subject to this invariant, institute i chooses the best remaining item and adds it into the bin (i, j), 
which has the maximum available space in the current phase. This unavoidably raises another issue: 
how can we be sure that there is at least one remaining item for institute i to put into the bin 
without violating Invariant CI We will address this issue in our proof. 

Theorem 24. Let x € Pfsm- Let Mij be the set of applicants assigned to bin {i,j) at the end of 
any step of the packing procedure and aij be the lowest-ranking applicant of institute i in bin {i,j) 
(implying Xiai^ is on top of bin {i,j)). Then 

(i) In any step, suppose that the algorithm is examining bin {i,j)- Then institute i can find at least 

one item in its remaining items to add into bin {i,j) without violating Invariant C; 
(a) For all bins {i,j), x{Mij\{aij}) + x{T{i,aij)) = x{Mij) + x{T{i,aij)\{{i,aij)}) = 1; 
(Hi) At the end of any step, institute i can organize a comb K{i, S{A^)) where Al is composed of appli- 

cants in {aij.ljfif)! so that x{K{i,SiA')) = ^if:?' ^(M^y ) + X:if=?' x(r(i, a^jOVK^, = 

\Ei{x)\; 

(iv) At the end of any step, an item Xia is not put into institute i 's bins if and only if there exists a 
class G a{C{t)) so that \{a,,,>}f:}f^ n n L^J = g+(C:,). 

(v) If Xia is packed and Xiia is not, then i' >-a i; 

(vi) At the end of any phase, the Uij in all bins are distinct. In particular, for any applicant a who 
is matched under x, there exists some bin {i,j) such that a = aij. 

Proof. We first assume that (ii) holds and prove (i) . Observe that (ii) implies that given any applicant 
a € Ei{x), its corresponding item Xia, if already put into a bin, must be on its top and fills it completely. 
Since {i,j) currently has available space, at least one applicant in Ei{x) is not in institute i's bins yet. 
We claim that there exists at least one remaining applicant in Ei{x) that can be added into bin (i, j). 
Suppose not. Let the set of applicants in Ei{x) that are not put into i's bins be G. Given any applicant 
a & G, there must exist some class C^ € a(C(i)) for which | Ui<j'<|_Bi(x)| j'j^j ^ ~ Q^i^l)- 
B be the set of classes C|. that contains at least one applicant in G and | Ui<j'<|_Ei(a;)| j'^j '^^fcl ~ 
q~^{Cl.). Let G' be {Ei{x)\G)\[J(ji^^i^^-jGl, the subset of applicants in Ei{x) that are already put 
into the bins but not belonging to any class in ^(B). Note that none of the applicants in G' can be 
in the bin (i, j). Thus, by counting the number of the bins minus (i, j), we have 

\Ei{x)\ 

\E,ix)\-l>\G'\+ Yl I U HrnCl\ = \G'\+ Yl l^iCl) 
C^e5R(B) i'=ij¥i Cl&R{B) 

Note that all applicants in Ei{x)\G' are in some class in (either they are already put into 

the bins or not). Then by the pigeonhole principle, there is at least one class C\ G K(-B) for which 
\{Ei{x)\G') n Ci\ > g+(C^), contradicting Lemma [22l 
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We now prove (ii)-(vi) by induction on the number of phases. In the beginning, (ii)(v)(vi) holds by 
Lemma [23r ii) (iii) . (iii)(iv) hold by setting ^* := Hi{x) and observation Definitionl21land Lemma l23( ii). 
Suppose that the theorem holds up to Phase t. Let a be the maximum available space in Phase t + 1. 
Suppose that the algorithm is examining bin (i, j) and institute i chooses item Xia to be put into this 
bin. From (vi) of the induction hypothesis, applicant a is on top of another bin where i' ^ i, 

in the beginning of phase t + 1. Then by (ii)(v) of the induction hypothesis, 

x{T{i, a)) < x{T{i', a)) - Xi>a = 1 - x{Mi'j') < a, (7) 

where the last inequality follows from our assumption that in Phase t + 1, the maximum available 
space is a. Note also that 

x{T(i,a)) = a, then G Lt+i (bin is also examined in Phase t+ 1). (8) 

Assume that A* is a tuple composed of applicants in ^' • For our induction step, let 

:= ]4*|"» Ja, the new set of items on top of i's bins after a is put on top of Uij. 

We first prove (iv). Since Xia is not put into the bin before this step, by (iv) of the induction 
hypothesis, there exists some class C*; G a{C{i)) for which D C\i n L!_„| = (?'''(C*J. Let C*; be 
the smallest such class. Since Xia is allowed to put on top of Xiai,j , ajj >-i a and aij G C*;, otherwise, 
Invariant C regarding q^iCl^j) is violated. 

Now we show that all other items xia' fulfill the condition stated in (iv). There are two cases. 

— Suppose that xia' is not put into the bins yet. 

• Suppose that Ojj- '^i a' >-i a. We claim that it cannot happen that for all classes C*,^ G a'(C(i)), 

n C*,j nL^^, I < q^{Cl^,^), otherwise, A*|"a' is still feasible, in which case institute i would 
have chosen Xia', instead of Xia to put into bin a contradiction. 

• Suppose that aij >-i a >-i a! . By (iv) of the induction hypothesis, there exists a class C*,^, G 
a'(C(i)) for which \K r\C\,^,r\\}^^, \ = g+(C*,,/). If C*/;/ (t C^i, it is easy to see that |A*nC*,,,n 
Ka'\ = q'-iCan')-^ if Ci^'i' C Ci^, then C^, G a'(C(z)) and we have \A^ n Q n hl^,\ = q+{Cl{). 
In both situations, the condition of (iv) regarding xia' is satisfied. 

— Suppose that xia' is already put into the bins. It is trivial if a' >-i a, so assume that a >-{ a'. We 
claim that none of the classes C*/^ G a'(C(i)) can be a subclass of (7*^ or itself. Otherwise, 

G a'{C{i)), and we have g+(C*;) = |X n C^, n L^] > \T n C*; n L!_^,|, a contradiction to 
(iv) of the induction hypothesis. Now since for every class C*,^ G a'(C(i)), we have C*,^ ^ C*;, we 

have \A^ D C*,^ n IL!_q,| = \A^ n C*,^ n L^, I < g+(C*,J, where the strict inequality is due to the 
induction hypothesis. 

We notice that the quantity Yl\^=i^^ x{Mij') is exactly the sum of the shaft x{S{A^)) (before Xia 
is added) or x{S{A^)) (after Xia is added) by observing (iv). Below let x{Mij) and x{Mij) denote the 
total size of the items in bin (z, j) before and after Xia is added into it. So x{Mij) = x{Mij) + Xia- 
Now we can derive the following: 
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\Ei{.x)\ 

x{Kii, S{A'))) = x{S{A')) + x{T{i, a)\{{i, a)}) + 

_ \Er{x)\ 

= x{Mij) +Xia + x{T{i,a)\{{i,a)}) + ^ x{Mij>) + x(r(i, aij/)\{(i, ajj/)}) 

= x{Mij) + x{T{i, a)) + \Ei{x)\ — 1 (by (ii) of the induction hypothesis) 
> \Ei{x)\ (by Constraint ([3])) 

For the above inequality to hold, 

x(Mij) + x(T{i,a))>l. (9) 

Since x{Mij) = 1 — a and x(T{i,a)) < a by Inequality ([7]), Inequality Q must hold with equality, 
implying that x{K{i, S{A'^))) = \Ei{x)\, giving us (iii). 

Since institute i puts Xia into bin the "new" Mij and the "new" a^j- (=a) satisfies 

x{Mij) + x{T{i,a)\{{i,a)}) = 1. 

This establishes (ii). (v) follows because Inequality ([7]) must hold with equality throughout. There- 
fore, there is no institute i" which ranks strictly between i and i' and Xj/z^ > 0. 

Finally for (vi), note that x{T{i, a)) = a if the item Xia is put into some bin in Phase t + l. All such 
items are the least preferred items in their respective "old" bins (immediately before Phase t + 1), it 
means the items on top of the newly-packed bins are still distinct. Moreover, from ([8]), if a bin 
is not examined in Phase t -\- 1, then its least preferred applicant cannot be packed in phase t -|- 1 
either. □ 

We define an assignment based on a number a G [0, 1) as follows. Assume that there is a line of 
height a "cutting through" all the bins horizontally. If an item Xia whose position in i's bins intersects 
a, applicant a is assigned to institute i. In the case this cutting line of height a intersects two items 
in the same bin, we choose the item occupying the higher position. More precisely: 

Given a G [0,1), for each institute i G X, we define an assignment as follows: /i"(i) = {a : 
1 - x{T{i, a)) < a < 1 - x{T{i, a)) + Xia}. 

Theorem 25. The polytope determined by Constraints |2])-([3) is integral. 

Proof. We generate uniformly at random a number a G [0, 1) and use it to define an assignment /i". 
To facilitate the discussion, we choose the largest a' < a so that /i" = /^". Intuitively, this can be 
regarded as lowering the cutting line from a to a' without modifying the assignment, and 1 — a' is 
exactly the maximum available space in the beginning of a certain phase I during the execution of our 
packing algorithm. Note that the assignment is then equivalent to giving those applicants (items) 
on top of institute z's bins to i at the end of phase I. 

We now argue that fj," is a stable matching. First, it is a matching by Theorem [2^vi). The 
matching respects the quota of all classes since Invariant C is maintained. What remains to be 
argued is the stability of fi". Suppose, for a contradiction, {i,a'^) is a blocking pair. We consider the 
possible cases. 
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— Suppose that x^^4, > and x^^,/, is not put into the bins yet at the end of Phase I. Then by 
Theorem I24(iv) and the definition of blocking pairs, {i,a'^) cannot block 

— Suppose that XjQ0 > and Xj^,/, is already put into the bins at the end of Phase /. If fi°'{a'^) = i, 
there is nothing to prove. So assume //"(a*^) 7^ i and this means that the item x^^^,/, is "buried" 
under some other item on top of some of i's bins at the end of Phase I. Then by Theorem I24lfv). 
a'^ is assigned to some other institute ranking higher than i, contradicting the assumption that 
{i,a'^) is a blocking pair. 

— Suppose that Xj^^ = 0. There are two subcases. 

• Suppose that for each of the classes C*^^ E a'^{C{i)), H C^^J < q^{C^^^^). Then we 
can form a new feasible tuple /x°(i)|a'^. It can be inferred from the definition of the shaft 
that x{S{^°^{i)\a'f')) < x{S{fi°'{i)). Moreover, by Theorem [2^iii) . we have x{K{i, S{fi"'{i))) = 
\Ei{x)\. Now by Constraint ([3]), 

\Ei{x)\ + l<x{K{i,S{fi''{i)\a'^))) 

<x(5(A."(i))+x(r(z,a*)\{(i,a*)})+ J] x{T{i,a)\{{i,a)}) 

aefJ.°' 

= x{K{i, S{^%i)))) + x(T(z, a<^)\{(z, a^))) 
= \Ei{x)\+x{T{i,a'i>)\{{i,a'*>)}). 

As a result, x(T(i, a'^)\{(i, a*^)}) = 1, implying that ^°'{a'^) l^^^ i, a contradiction to the 
assumption that {i,a) blocks /x". 

• Suppose that there exists a class C*^^^ G a^{C{i)) for which n C^.^J = Let 
C^^^^ be the smallest such class. By definition of blocking pairs, there must exist an applicant 

G ^a'l'i'P ranks lower than a'^'. Choose to be the lowest ranking such applicant 

in We make the following critical observation: 

x{S{^,^{^fa^))<x{S{^^%^))-x,,,. (10) 

To see this, we first argue that given an item Xia > 0, if it does not contribute to the shaft 
S{iJ,'^{i)), then it cannot contribute to shaft S{fj,"{i)\"''' a'^) either. It is trivial if a )~i a) . So 
assume that >-i a. First suppose that a C'*^;^- Then given any class C*j G a(C(i)), 
n Cit n hiJ = |;u"(i)|"^a'^ n Ci^ n L*^„|, and Theorem [Mtiv) states that there is a class 
C*; G a{C{i)) such that n n L^l = g+(C*;). Secondly suppose that a G C^^^^^. 

Observe that = \^Ji''{i)f a't> n C[^^, D l^lJ = \fi'f'{i)\^' D C\,^, H L^] (the first 

equality follows from the choice of o^). In both cases, we conclude that Xia cannot contribute 
to the shaft S'(/i'^(i)|"^a'^). The term x^^t does not contribute to the shaft S{fi'f'{i)\"''' a'^') by the 
same argument. Now using Constraint ([3]), Theorem IMT iii) . and Inequality (fTO]) . we have 

\m{x)\<x{K{i,S{fi"{i)fa''^))) 

< x(5(/.-(i))) - x,,t + xiTii, a<^)\{(., a<^)}) + ^ x(T(., a)\{(., a)})) 

ae^l'^{i)\{a^ 

= \Ei{x)\ - x{T{i,a^)) + x{T{i,a't')). (Note that x^^,^ = 0). 

Therefore, 
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x{T{i, a'^)) > x{T{i, a^)) > 1 - a' > 1 - a. 

So f.i°'{a'^) )^^4, i, again a contradiction to the assumption that {i^a"^) blocks 
So we have established that the generated assignment /i" is a stable matching. Now the re- 
maining proof is the same as in [23]. Assume that ;u"(z,a) = 1 if and only if applicant a is 
assigned to institute i under /i". Then 

Exp[iJ°'{i,a)] = Xia- 

Then Xia = Jq fJ-'^ih a)da and x can be written as a convex combination of /x" as a varies over 
the interval [0, 1). The integrality of the polytope thus follows. □ 

4.1 Optimal Stable Matching 

Since our polytope is integral, we can write suitable objective functions to target for various optimal 
stable matchings using Ellipsoid algorithm [10]. As the proposed LP has an exponential number of 
constraints, we also design a separation oracle to get a polynomial time algorithm. The basic idea of 
our oracle is based on dynamic programming. 

4.2 Median-Choice Stable Matching 

An application of our polyhedral result is the following. 

Theorem 26. Suppose that in the given instance, all classifications are laminar families and there 
is no lower bound, q~{Cj) = for any class Cj. Let ni, ^2, ■ ■ ■ , l''k be stable matchings. If we assign 
every applicant to his median choice among all the k matchings, the outcome is a stable matching. 

Proof. Let x^^ be the solution based on fit for any 1 < t < k and apply our packing algorithm on the 

fractional solution x = ^* . Then let a = 0.5 and lu,^'^ be the stable matching resulted from the 

cutting line of height a = 0.5. We make the following observation based on Theorem 1241 

Suppose that applicant a is matched under x and those institutes with which he is matched 
are ii, i2, • • • , ik', ordered based on their rankings on a's preference list. Assume that he is 
matched to it rit times among the k given stable matchings. At the termination of the packing 
algorithm, each of the items Xj^a, 1 < Z < A;', appears in institute i/'s bins and its position is 
from Etl f to Et=i f • 

Now /i"'^ gives every applicant his median choice follows easily from the above observation. □ 

Using similar ideas, we can show that an applicant-optimal stable matching must be institute- 
(lexicographical)-pessimal and similarly an applicant-pessimal stable matching must be institute- 
(lexicographical)-optimal: by taking x as the average of all stable matchings and consider the two 
matching fi'' and /i^""^ with arbitrary small e > 0. Hence, it is tempting to conjecture that the 
median choice stable matching is also a lexicographical median outcome for the institutes. Somehow 
surprisingly, it turns out not to be the case and a counter-example can be found in the appendix. 
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4.3 Polytope for Many-to-Many "Unclassified" Stable Matching 



In the many-to-many stable matching problem, each entity e I L) A has a quota Q{e) S Z+ and 
a preference over a subset of the other side. A matching fj, is feasible if given any entity e (z I U A, 
(1) |^(e)| < (5(e), and (2) ^(e) is a subset of the entities on e's preference list. A feasible matching 
H is stable if there is no blocking pair (i,a), which means that i prefers a to one of the assignments 
or if < Q{i) and a ^ and similarly a prefers i to one of his assignments /x(a), or if 

|//(a)| < Q{a) and i ^ fi{a). 

We now transform the problem into (many-to-one) LCSM. For each applicant a G we create 
Q{a) copies, each of which retains the original preference of a. All institutes replace the applicants 
by their clones on their lists. To break ties, all institutes rank the clones of the same applicant in an 
arbitrary but fixed manner. Finally, each institute treats the clones of the same applicant as a class 
with upper bound 1. It can be shown that the stable matchings in the original instance and in the 
transformed LCSM instance have a one-one correspondence. Thus, we can use Constraints ([I])-(l5]) 
to describe the formeiH. 

5 Conclusion and Future Work 

In this paper, we introduce classified stable matching and present a dichotomy theorem to draw 
a line between its polynomial solvability and NP-completeness. We also study the problem using the 
polyhedral approach and propose polynomial time algorithms to obtain various optimal matchings. 

We choose the terms "institutes" and "applicants" in our problem definition, instead of the more 
conventional hospitals and residents, for a reason. We are aware that in real-world academics, many 
departments not only have ranking over their job candidates but also classify them based on their 
research areas. When they make their hiring decision, they have to take the quota of the classes into 
consideration. And in fact, we were originally motivated by this common practice. 

CLASSIFIED STABLE MATCHING has happened in real world. In a hospitals/residents matching 
program in Scotland, certain hospitals declared that they did not want more than one female physician. 
Roth [16] proposed an algorithm to show that stable matchings always exist. 

There are quite a few questions that remain open. The obvious one would be to write an LP 
to describe LCSM with both upper bounds and lower bounds. Even though we can obtain various 
optimal stable matchings, the Ellipsoid algorithm can be inefficient. It would be nicer to have fast 
combinatorial algorithms. The rotation structure of Gusfield and Irving [11] seems the way to go. 

Acknowledgments 

I thank Peter Winkler for suggesting to me the idea of classified stable matching and Alvin 
Roth and Jay Sethuraman for several fruitful discussions. The idea of reducing many-to-many stable 
matching to LCSM was suggested by Jay Sethuraman. The comments of SODA 10 reviewers have 
been extremely helpful. 

^ This cloning technique works for the many-to-one hospitals/residents problem, i,.e, we can clone the hospitals and 
use the original one-one stable marriage polytope [24] to describe that of the former [20l22j . However, this trick does 
not work here for many-to-many matching. The reason is that the same applicant cannot be assigned to an institute 
multiple times. And that is why our classifications help to resolve this problem. 



23 



References 

1. Abeledo, H. G., and Rothblum, U. G. Stable matchings and linear inequalities. Discrete Applied Mathematics 
54(1) (1994), 1-27. 

2. Abraham, D., Irving, R., and Manlove, D. Two algorithms for the student-project allocation problem. Journal 
of Discrete Algorithms 5(1) (2007), 73-90. 

3. Ba'iou, M., and Balinski, M. The stable admissions polytope. Mathematical Programming, series A 87 (2000), 
427-439. 

4. Biro, P., Fleiner, T., Irving, R., and Manlove, D. The college admissions problem with lower and common 
quotas. University of Glasgow, Computing Science Department Research Report, TR-2009-303, 2009. 

5. Fleiner, T. Some results on stable matchings and fixed points, 2002. TR-2002-08, EGRES, December. 

6. Fleiner, T. A fixed-point approach to stable matchings and some applicants. Mathematical Operations Research 
28(1) (2003), 103-126. 

7. Gale, D., and Shapley, L. College admissions and the stability of marriage. American Mathematical Monthly 
69(1) (1962), 9-15. 

8. Gale, D., and Sotomayor, M. Some remarks on the stable matching problem. Discrete Applied Mathematics 11 
(1985), 223-232. 

9. Garey, M., and Johnson, D. Computers and Intractablility. Freeman, 1979. 

10. Grotschel, M., Lovasz, L., and Schrijver, A. The ellipsoid method and its consequences in combinatorial 
optimization. Combinatorica 1(2) (1981), 70-89. 

11. GuSFlELD, D., AND Irving, R. The Stable Marriage Problem. The MIT Press, 1989. 

12. Hamada, K., Miyazaki, S., and Iwama, K. The hospitals/residents problem with quota lower bounds. In Match- 
Up: Matching Under Preferences-Algorithms and Complexity, Satellite workshop of ICALP 2008 (2008), pp. 55-66. 

13. Klaus, B., and Klijn, F. Median stable matching for college admissions. International Journal of Game Theory 
34(1) (2006), 1-11. 

14. Knuth, D. Mariages stables et lews relations avec d'autre problemes combinatoires. Les Presses de I'universite de 
Montreal, 1976. 

15. Roth, A. On the allocation of residents to rural hospitals: a general property of two-sided matching markets. 
Econometrica 54 (1986), 425-427. 

16. Roth, A. A natural experiment in the organization of entry level labor markets; Regional markets for new physicians 
and surgeons in the U.K. American Economic Review 81 (1991), 415-440. 

17. Roth, A., Rothblum, U. G., and Vate, J. H. V. Stable matchings, optimal assignments and linear programming. 
Mathematics of operation research 18 (1993), 808-828. 

18. Roth, A., and Sotomayor, M. Two-sided matching: A study in game-theorectic modeling and analysis. Cambridge 
University Press, 1990. 

19. Rothblum, U. G. Characterization of stable matchings as exterme points of a polytope. Mathematical Programming 
54 (1992), 1475-1480. 

20. Sethuraman, J. Private communication, 2009. 

21. Sethuraman, J., and Teo, C.-P. A polynomial-time algorithm for the bistable roommates problem. Journal of 
Computer and System Sciences 63(3) (2001), 486-497. 

22. Sethuraman, J., Teo, C.-P., and Qian, L. Many-to-one stable matching: Geometry and fairness. Mathematics 
of Operations Research 31(3) (2006), 581-596. 

23. Teo, C.-P., AND Sethuraman, J. The geometry of fractional stable matchings and its applications. Mathematics 
of Operations Research 23(4) (2001), 874-891. 

24. Vate, J. H. V. Linear programming brings martial bliss. Operation Research Letters 8 (1989), 147-153. 



24 



A An Example for Section 12.21 



In contrast to the generalized rural hospitals theorem in LCSM, if some institutes use intersecting 
classes, stable matching sizes may differ. Figure [3] is an example. 



Institute Preferences Classifications Quota 

iiMiazaa Cl = {a^, az}, Ck ^ {ai, a-^} Q{ii) = 2, q+ (Cl) = 1, q+ (Ck) = 1 

h-.a-iaxa^a^ Cl = {aa.ai}, C| = {aa.ag}, C| = {a^.a^) Q(i2) = 2, g+(C?) = 1, <?+(C|) = 1, q+(C|) 



Applicant Preferences 

ai:iii2 

a2:iii2 

03:1112 

04:12 



Stable Matchings 

HA = {(ii;ai), (12; 02)} 

Ms = {(ii; 02,03), (12; oi, 04)} 

Fig. 3. An example of stable matchings of different sizes. 



B Missing Proofs of Section [3] 

In this section, we prove Theorem [161 We assume that the set of posets P = {Pi, P2,- " ,Pk} contains 
a poset which is not a downward forest. Moreover, we assume that there is no lower bound on the 
classes. 

Without loss of generality, we assume that Pi is not a downward forest. Such a poset must have 
a "V." By definition, there exists institute i whose class inclusion poset P{i) is isomorphic to -Pi. 
This implies that institute i must have two intersecting classes in C(i). In the following, we will 
present a reduction in which all institutes use at most two classes (that can be intersecting). It is 
straightforward to use some dummy institutes and applicants to "pad" our reduction so that every 
poset Pj G P is isomorphic to some class inclusion poset of the institutes in the derived instance. Our 
reduction is from ONE-IN-THREE-SAT. We will use an instance in which there is no negative literal. 
(NP-completeness still holds under this restriction [9].) 

The overall goal is to design a reduction so that the derived P-CLASSIFIED STABLE MATCHING 
instance allows a stable matching if and only if the given instance = ci A C2 A • • • A Cfc is satisfiable. 
We will build a set of clause gadgets to represent each clause Cj. For every pair of literals which belong 
to the same clause, we create a literal-pair gadget. Such a gadget will guarantee that at most one 
literal it represents can be "activated" (set to TRUE). The clause gadget interacts with the literal- 
pair gadgets in such a way that if the clause is to be satisfied, exactly one literal it contains can be 
activated. 

Literal-Pair Gadget Suppose that and x^, both belong to the same clause Cj. We create a gadget 
T--, composed of four applicants {a^^}^^^ U {a^/ ^1^=1 and two institutes whose preferences 

and classifications are summarized below. 
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H- 4,2 >- <i >- 4,2 >- 4,1 >^ ^iH) c'l = {4i'42}> ci' = {4i'4,i} 



4: all ^42 ^4,1 ^4,' 



g(7/) = 2, q+{C^) = l,q+{C^) 

Ci = {ail.ao} 
Q(/^,) = 2, g+(c'0 = l 



We postpone the explanation of the F and ^ functions for the time being. We first make the 
fohowing claim. 

Claim B: Suppose that in a stable matching the only possible assignments for {4 !> 4 2' 4' i' 4' 2} 
are {Ij , /^,}. Then there can only be three possible outcomes in fj,. 

1. /i(4,i) = ^iJ /"(4.2) = ^i'^ A*(4',i) = ^i'' /^(4'.2) = ^i - ^'^^^ '"^^ ^« activated 
while x^/ remains deactivated.) 

2. /i(4,i) = ^i', /^(4,2) = 4' A*(4',i) = 4' /^(4,2) = 4" (^^ ^^^^ ^^s^' '^^ ^i' activated 
while remains deactivated.) 

3. Ai(4,i) = 4" Ai(4,2) = 4' /"(4',i) = 4" /^(4',2) = 4- (I'^ t^^^^ ^^^^e, we say both Xi and Xj/ 
remain deactivated.) 

Claim B can be easily verified. Note that the case fJ-iaj^) = If, nial^) = 4" /^(4' i) ~ 4' 

^(4' 2) ~ 4' ^^^^ '^'-'^ happen due to the quota q~^{C2^ )■ This case corresponds to the situation that 
Xi and Xi' are both activated and is what we want to avoid. 

We now explain how to realize the supposition in Claim B about the fixed potential assignments 
for {44?=i {4' *}*=! ^ stable matching. It can be easily checked that if 4i is matched to some 
institute in 1^(4 1), or either of {4 1542) unmatched; or if either of {aj, ^,aj, 2} is unmatched, 
then there must exist a blocking group involving a subset of {4 > 4" i4 t}?=i44' *}*=! }• However, 
the following matching //"^ can happen in which 4' 1 is matched to some institute in i~'(4' 1) but there 

is no blocking group : /i'^(4,i) = 4, /U^(4,2) = f^H4,2) = 4' /^'^(4m) ^ ^(4,i)Ill 

To prevent the above scenario from happening (i.e., we want /i*^ to be unstable), we introduce 
another gadget T^, associated with Ij, to guarantee a blocking group will appear. We now list the 
preferences and classifications of the members of below. 



,4: Ii,4 >- Ii,l >- I', 

,5 



1,3 Ii,4 I', 
,3- ^lA >~ Ii,Z ^ I', 



■ Ii,2 >- Ii,4 >- I', 



,3^4 
,2^4 

,1^4 
,2^4 

,3^4, 



i,2 ^ Ii,4 >~ Ii,3 ^ Ii,l 



73 , . . ■ 

{41.42.43}' 
1 



-li,2- "i,4 



>- aiji >- <,2 >- "Is >- all >- "Is C'l 

Q(Ji2) = 2, q+(C7^''^) = l, 

>- 45 ^ "4 >- o-is >- «li ^ "I2 C'l''"" = {«4.^l2.al3}' ^'2'''' = {al3.al4.al5} 



g(7l4) = 2, = 1, g+(C2-^) = 1 





al4. 


n 


(c4^) 


= 1 




= {"Is. 


al4. 










= 1 




= {"Is. 


al4. 




7^ 

(^2-4 


= 1 





It can be verified that if ^ is matched to some institute in r{a'^, j), the above assignment is the only possibihty 
that no blocking group arises. 
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The above instance has the following features, every one of which is crucial in our construction. 

Feature A does not allow a stable matching; more importantly, if in the given matching /x, 
< 2, then there exists at least one blocking group which is not of the form {li^^al ^,al y), 
where C {a],J?=i- 

Feature B When the institute i is "removed" from the instance (i.e. i is struck from the 
preferences of the applicants {a^^}^^^), there exists at least one stable matching. For instance, 
the following matching /j, is stable, if lli is removed: /i(/i2) — (^le'^is)' /^(-^la) ~ (^i5'^i2)' 

Feature C Each of the three institutes {Iit}t=2 ^^^s exactly two intersecting classes; more impor- 
tantly, institute ll i does not use any classification at all. 

Our idea is to let the institute in the gadget T--, "play the role" of I^ i in the gadget of and 

add the other members in into gadget T--,. To be precise, let 1^(1/ ) = aj^ ^ 

cij^ y all (^^^ same preference list of l]^) and add the other members of into the gadget T^-, 
without modifying their preference lists (and classifications). We now explain how the above features 
of tI help to realize the supposition in Claim B about the potential assignments of the applicants 
{a\^,a\i in a stable matching. 

1. In a matching , suppose that institute is only assigned a^^ while a^, ^ is assigned to some 
institutes in r{a'-, ^ (the problematic case we discussed above). As a result, institute if can 
take one more applicant from the set {a^^jj^j^. By Feature A, there must exist a blocking group 
involving the members in . More importantly, this blocking group need not be composed of If 
and two applicants from {a]^}^^]^. 

2. In a matching /i*^, suppose that institutes if is assigned two applicants from the set {al^, a^, j}f=i- 
Then l\ ^ can be regarded as being removed from the instance T\ . And there exists a stable 
matching among the other members of the instance T\. This explains the necessity of Feature B. 

3. Finally, since If already uses two intersecting classes, ^ should not use any more classes. This 
explains the reason why Feature C is necessary. 

We have left the functions r[ali) and r{ali ^) unexplained so far. They contain institutes belonging 
to the clauses gadgets, which will be the final component in our construction. 

Clause Gadget Suppose that Cj = x{ V V Xg. We create a clause gadget Tj composed of two 
institutes {/j}^^^ and six applicants {al}^^^. Their preferences and classifications are summarized 
below. 

We now explain how the A functions in the clause gadgets interact with the F functions in the 
literal-pair gadgets. The former is composed of applicants in the literal-pair gadgets while the latter is 
composed of institutes in the clause gadgets. Our intuition is that the only possible stable matchings 
in the clause gadgets will enforce exactly one of its three literals to be activated. To be precise, let 
7r(X) denote an arbitrary order among the elements in the set X. Then: 
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ai: i{ y ii 
ai:Uyii 
ai: ifyU 



!{: a^y a{y a^y A{x{) y a^ y A{xi) y a^^y A{xi^) y a\ 
= {a{,ai,,ai,al,A{x{)}, C^^ = {ai,ai,al, A{xi)} 



al y ai y a{ y al y ai y tig 



Ci^ = {a{,ai,al,al}, C^" = {aj[,a^,a^}. 
Q{Ii) = 3, q+{c'j) = 2, = 2 



A{x{) = 7r({a{j if a(\ G T^'^, and xf = x{}), for t e {1, 2, 3}. 
-^('^i i) = '''"({-^1 I if clause Cj/ = x\ V V contains the literal x^.}) 

Claim C Suppose there exists a stable matching /z in the derived instance. Then there can be 
only three possible outcomes for the members in the clause gadget Tj. 

1. = ('^5) ^i) ^^3)1 /"(-^2) ~ ('^ei' '^2' ^i)- Moreover, in this case, X2 must be activated while 
x\ and X3 must remain deactivated. 

2. = (dg, 02)^4), ^{iD ~ (^6''^!' ^5)- Moreover, in this case, Xg must be activated while 
x\ and X2 must remain deactivated. 

3. = (d2; 05,03), ^i{I2) = (01,05,04). Moreover, in this case, x\ must be activated while 
X2 and X3 must remain deactivated. 

Proof of Claim C. We have argued previously that in a stable matching ji, the applicants in the 
literal-pair gadgets cannot be matched to the institutes in the clause gadget. So both and can 
only have applicants from the set {d^jj^i in /i. We first make the following observation. 

In a stable matching ^, one of the applicants in the following three pairs {01,02}, {03,04}, 
{d5,dg} must be assigned to !{; while the other applicant in the above three pairs must be 
assigned to Ig. 

The above observation can be easily verified. So there are eight possible combinations in Except 
the three combinations listed in the claim, the other five will cause blocking groups to arise. 

1. = (4,^14)^ t^i^i) = Then {ll; ^t{l()\^2al) blocks fi. 

2. //(/i) = (d{,d;^,d{), fJ-{ii) = {ai,di,ai). Then again fi{ii)\°'2ai) blocks /j,. 

3. /i(ii) = (05,02,03). Note that the tuple //(I^) = {4^'^i^4) feasible. So /^(/j) = (dg,dj). 
But then {I2,n{l2)\4) blocks fi. 

4. ^(/i) = (05,01,04), ^(12) = (og, 02,03). In this case, institute can replace 04 with an applicant 
in vl(x2) or an applicant in vl(x3) to have a blocking group, unless both Xg and Xg are activated. 
However, this is impossible, since we have a literal-pair gadget T23 and by Claim B, at most one 
of X2 and X3 can be activated. So we have a contradiction. 

5. /-f(/i) = (02,05,04), ^{12) = (01,05,03). Then the existence of a blocking group involving !( and 
some applicants in yl(xi) Uyl(x3) is guaranteed. The reason is similar to the previous case: xj and 
Xn cannot be both activated. 
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Finally, we remark that the three possible outcomes in listed in the lemma will guarantee that 
exactly one of the three literals in clause cj can be activated. The reason is again the same as in the 
last two cases that we just explained. This completes the proof of Claim C. □ 

Now by Claim C, we establish Theorem [TBI 

C Missing Proofs of Section [4] 

Lemma 17. In LCSM, if there is no lower hound, i.e., given any class C'j, q~{Cj) = 0, then a stable 
matching as defined in Definition\^can he equivalently defined as follows. A feasible matching /i is sta- 
ble if and only if there is no blocking pair. A pair (i, a) is blocking, given that fi{i) = {aii,ai2, • • • , aik), 
k < Q{i), if 

- i >-a 

- for any class C^t £ a{C{i)), n n C^J < q+iCit). 

Proof. If we have a blocking group {i; g), institute i and the highest ranking applicant in g\fJ-{i) must 
be a blocking pair. Conversely, given a blocking pair (z;a), assuming that |/u(z)| = Q{i) (the case 
that |/u(i)| < Q{i) follows a similar argument), we can form a blocking group (i; /i(i)|"^a), where 
is chosen as follows: (1) if there exists a class C*j S a(C(i)) such that n C*^| = ^^(C*^), choose 
the smallest such class C*^ G a{C{i)) and let be the lowest ranking applicant in n C*^; (2) 
otherwise, is simply the lowest ranking applicant in fi{i). □ 

Lemma 19. Every stable matching solution x satisfies the comb inequality for any comb K{i, S{A'^)): 
x{K{i,S{A')) = + ^iTii,aj)\{i,aj}) > \A'\. 

We use the following notation to facilitate the proof. Give a tuple A', we define yia as follows: 
_ J 1 either a £ A\x{T{i, a)) = 1; or a A\ Xia = 1, and (i, a) £ S{A^); 

Via — ^ n 

1^ o.w. 

Let y(q) = E 

aeL'nc^ This quantity indicates how much a class Cj contributes to the comb 
value x{K{i, S{A'^))). Thus, if [/ is a set of classes in C{i) partitioning L*, then x{K{i, S{A'^))) = 

Proof. We prove by showing that if x{K{i, S{A'^))) < \A^\, there exists a blocking pair {i,a^), where 
G A^. We proceed by contradiction. First note that there exists a non-empty subset G C of 
applicants a for whom x{T{i,a)) = 0, otherwise, x{K{i, S{A^))) > \A^\, an immediate contradiction. 
For each applicant a £ G, there must exist a class C*^ £ a{C{i)) for which Yla'^L^ nc\ ~ 1~^(^ai)^ 
otherwise, (i, a) is a blocking pair and we are done. Now for each applicant a £ G, choose the smallest 
class C*; for which X^^j/g^i f^Qi^ Xia' = q~^iGl^i) and denote this class as Ga- We introduce a procedure 
to organize a set U of disjoint classes. 

Let G be composed of ai, a2, • • • , a|c| ordered based on their decreasing rankings on L* 
For i = 1 To \G\ 

if ai £ C £ U, then do nothing 

else U := U\{G\G £U,G C Ca^} //Gai may be a superclass of some classes in U 
U ■.= UU {CaJ. // adding Ca^ into U 
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Claim The output U from the above procedure comprises of a disjoint set of classes containing 
all applicants in G, and for each class G U, y{Cj) > q^{Cj). 

We will prove the claim shortly. Now 

xiK{i,s{A'))) = + l^^{uc^6^q}| > i^icj) + l^^{u^.ec/q}l > m 

C](iU ' CjGf/ 

a contradiction. □ 
Proof of the Claim. It is easy to see that the classes in U are disjoint and contain all applicants 

in G. Below we show that during the execution of the procedure, if Cj € U , then ?/(Cj) > q+(Cj). 

We proceed by induction on the number of times U is updated. In the base case U is an empty set 

so there is nothing to prove. 

For the induction step, assume that a; is being examined and C^j is about to be added into U . 

Observe that even though X^^jg^i ^jj^ Xia = q~^{Cai), there is no guarantee that if Xia = 1, then 

yia = 1 for each a S nCa, . The reason is that there may exist some class Cj G a(C(i)) for which 
\A' n n hlJ = g+(Cj) and a A\ Then {i, a) is not part of the shaft x{S{A')) and y^a = 0. 

To deal with the above situation, we need to do some case analysis. Let B be the set of subclasses 
q of Ca, for which |A^nCjnL^J = Choose D to be the subclasses of C^, so that ^{BUU)UD 

partitions Ca, . We make three observations below. 

(i) for each class q G ^{B U U) and q G U, y(q) > g+(q) > Eae^^ncj ^ia- 

(ii) for each class q G D, if a G L!_^^ fl q and Xia = 1, then y^a = I. 

(iii) for each class q G U U) and q ^ ?7, then for each apphcant a G L^^^j^ n q n 
either a £ G and o G C G J7, or that a ^ G (implying that x{T(i,a)) = 1). Moreover, 

y(q) > SaeLt„ nc^ ^jfi 

(i) is because of the induction hypothesis and the feasibility assumption of x. (ii) follows from 
the fact that a ranks higher than a; and the way we define a class in D. For (iii), first notice that 
if q G ^{B U U) and q U, then such a class q must be part of ^{B) and q may contain 
some classes in U. Now suppose that aj G G n L'J_^^ but does not belong to any class in U. Then our 
procedure would have added the class Gq- into U before examining a;, a contradiction. To see the 
last statement of (iii), let G' be set of applicants in l^yai H Gj fl A* who do not belong to any classes 
in U. Then 



y(q)> E yiCi) + \G'\> J2 <l^iCi) + \G'\>q+{C^j)> Yl 



Xin 



where the first inequality follows from the first part of (iii), the second inequality the induction 
hypothesis, the third the fact that G] G ?R.{B) (thus |L'^^^ n G] n = g+(Gj)), and the fourth the 
feasibility assumption of x. 

Now combining all the three observations, we conclude that 

y(Ca,)= Yl 2/(q)+Ey(^j)^ E E ^- = ^^(q)' 

and the induction step is completed. □ 
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Lemma 20. Every stable matching solution x satisfies the following inequality for any class-tuple t*-; 



^ x{T{i,aij)\{i,aij})> ^ Xia (*) 

j 



Proof. We prove by contradiction. Suppose that in a given class-tuple t*- (*) does not hold. We will 
show that we can find a blocking pair (i,o^), where € t*-. Let the set of applicants a € with 
x{T{i,a)) = be G, a = J2a'eV nC' ^ ^' ^^^'^ ~ Yla'et^ ^ia'- assumption, at most a — 1 

J 

applicants a G have x(T(i, a)\{(i, a)}) = 1. Thus, 

|G| >g+(q) -/?-« + L (11) 



Claim: At least one applicant € G belongs to a sequence of classes CV^^ G a^(C(i)) such 
that if C^^^ C q, then Y.a'eV ,ncv ^m' < e+CC^t^). 

We will prove the claim shortly. Observe that given any class D C*, X^a'eL* nC' ^ ^'''(^fc)- 

as a > 0, some applicant a'^ € ranking lower than has Xj^.;, = 1 and Constraint ([2]) enforces that 
Sa'eL'nc^ — '/'''(C'fc)- Combining the above facts, we conclude that (i,a^) is a blocking pair. □ 

Proof of the Claim. We prove by contradiction. Suppose that for every applicant a € G, there 
exists some class C*^ G a{C{i)), Cl^ C C], and E„,gL!.„nc\ " l^i^at)- 

Let be the set of classes C C* such that contains an applicant a G G and Yla'ehi_ nc| ~ 
g"'"(C^) (which then will equal J2a'GV nC' "-^^^ *° Constraint ([2])). For each class G ^{B), 



x,a = q+{Cl)>\t)ncl\= Yl x,a' + \Gncll (12) 

where the first inequality follows from the definition of the class-tuple. Now we have 

q+{q)-a-P> Y E ^-'^ E \GnCl\ = \G\>q+iq)-a-P + l, 

a contradiction. Note that the first inequality follows from Constraint ([2]), the second inequality 
from (|12p . the equality right after is because every applicant in G belongs to some class in B, and 
the last inequality is due to (fTTj) . □ 

D Separation Oracle in Section 14.11 

It is clear that Constraints ([H) ([2]) (jj]) can be separated in polynomial time. So we assume that x 
satisfies these constraints and focus on finding a violated Constraint ([3]) and/or Constraint ([5]). 
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Separating Constraint ([3]) We first make an observation. For each institute i, it suffices to check 
whether ah the combs with exactly Q{i) teeth satisfy Constraint ([3]). To see this, suppose that there 
is a feasible tuple with less than Q{i) applicants and x{K{i, S{A^))) < \A^\. Then we can add 
suitable applicants into A to get a feasible tuple A^ with exactly Q{i) applicants. Noticing that 
x{S{A')) < x{S(A')), we have 

x{K{i, S{A^))) < x{S(A')) + J2 a)\{(^, a)}) + ^ x{T{i, a)\{{i, a)}) 

<\X\+ x{T{i,a)\{{i,a)}) 

< \A^\ + l^^l - \A'\ 
= \A% 

where the last inequality follows from our assumption that x satisfies Constraint ([T]). 

To illustrate our idea, we first explain how to deal with the case that the original classification 
C{i) is just a partition over L* (before we add the pseudo root class C|). We want to find out the tuple 

of length Q{i), whose lowest ranking applicant is , which gives the smallest x{K{i, If we 

have this information for all possible a^^, we are done. Note that because of our previous discussion, 
if there is no feasible tuple of length Q{i) whose lowest ranking applicant is , we can ignore those 
cases. 

Our main idea is to decompose the value of x{K{i, S{A^))) based on the classes and use dynamic 
programming to find out the combinations of the tooth-applicants that give the smallest comb values. 
More precisely. 

Definition 27. Assume that A^ C Cj, < | < (7+(Cj), and all applicants in Aj rank higher than 
a^. Let 

xiAi,a^)= J2 Xia+Y,xiTii,a)\{ii,a)}) 

a£V ^nCl{i,a)GS{A)) aGA) 

Z{Cj, Sj, a^) = minj^yj^i^(^Qi^^j^i^^g.x{Aj, a^). 

Note that this definition requires that if Xia contributes to x{A'j,a^), then a has to rank higher 

than , belongs to Cj, and the (i,a) is part of the shaft S{A'j). 

Suppose that we have properly stored all the possible values of Z{Cj, Sj,a^) and assume that 
G Cj,. Then for each class Cj ^ Cy, assume that < Sj < q~^{Cj) and for class Cj/, < sj' < 

q'^{Cj,) — 1, then the tuple A* whose lowest ranking applicant is a^, that gives the smallest comb 

value is the following one: 

The above quantity can be calculated using standard dynamic programming technique. So the 
question boils down to how to calculate Z{Cj, sj, a^). There are two cases. 
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— Suppose that Sj < q~^{Cj). Observe that all the positive Xia, where a ^ C"^- n L*^^| will contribute 
to the value x{A^pa)). So, to calculate Z{Cj,Sj,a^), we only need to find out the sj applicants a 
in Cj n with the smallest values x(T{i,a)\{{i,a)}). This can be easily done in polynomial 
time. 

— Suppose that Sj = q~^{C'!j). Different from the previous case, some of the positive Xia, where 
a G Cj n ILJ^^t will not contribute to x{A'j,a^). Our idea is to "pin down" the lowest ranking 
applicant E Aj. After we pin down a^, we know that the positive XiaS that contribute to 
x{A^j,a^) are those with a G n LJ^^j, while those with a G C* n n do not. So what 
remains to be done is to find out the sj — 1 applicants a in Cj fl ^l_^x with the smallest values 
x{T{i,a)\{{i,a)}). We then enumerate all possible G Cj n L*^^| and find out the applicant 
that gives the smallest -^(Cj, Sj, a^). The whole process can be done in polynomial time. 

We now explain how to generalize to the case that C{i) is a laminar family. In the previous 
simplified case that C(i) is a partition, we critically collect the following information: 

Suppose that Aj C Cj and < \Aj\ < q~^{Cj) and all applicants in Aj rank at least as 
high as and strictly higher than a"^. What is the choice of Aj so that x{S{A'j),a'^) + 
X^aeA' ^(^(^) '*)\{(^) '^)}) is minimized? 

The above question motivates the following definition. 

Definition 28. Assume that Aj C Cj, < |Aj| < q~^{Cj) and all applicants in Aj rank at least as 
high as a^, who in turn, ranks higher than a^. Furthermore, let x(Aj,a^) inherit the definition as 
defined in Definition\2'T\ We define 

^{^ji^ji^ 1^ ) ''^^''^Ay.A'jCC^.,\A'.\=Sj,all applicants in Aj rank at least as high as a^."^^^-?'^ ^' 

For the case that Z{Cj,Sj,a^,a'^) is not well-defined, e.g., there is no feasible tuple of length sj 
so that all applicants in Aj C Cj rank at least as high as a^, or a^, let Z{Cj, Sj,a^,a^) be an 

arbitrary large value. 

In the following, we also associate each Z{Cj, Sj, a^,a^) with the corresponding feasible tuple Aj that 
realizes this valu^. This helps us to identify the tuple that is violated in Constraint dSD- 

Theorem 29. For all classes Cj G C{i), we can correctly calculate all possible Z{Cj,Sj,a^,a^) in 
polynomial time, for all and combinations. 

Proof. We first remark that in the calculation, we indeed only need to find out all the values of 
Z{Cj,Sj,a^,a'^) for G Cj. For the case of ^ Cj, simply copy the value of Z{Cj,Sj,a^,a^) where 

a-t G Cj and is the lowest ranking applicant that ranks higher than a^. If there is no such applicant, 
let Z{Cj, Sj, a^, a^) be an arbitrary large value. 

We now give an inductive proof based on the height of class Cj in the tree structure of C{i). 
The base case is when Cj is a leaf class. Then the calculation of all Z{Cj,Sj,a^,a^) can be done in 
essentially the same way as we have shown in the case that C(i) is a partition. 

* If Z{Cj, Sj, a*, a^) is some arbitrary large value, let its corresponding tuple be (f> and indeed this will not matter. 
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For the induction step, let Cj be a non-leaf class and assume that a* e Cl G ciC^j). To calculate 
Z{Cj, Sj,a^,a^), we need to find out a feasible tuple Aj of size Sj, all of whose applicants rank at least 
as high as so that x{A'j, a^) is minimized. 

Observe that a feasible tuple Aj can be decomposed into a set of tuples Aj = Uc*gc(C') ^fe' '^here 
Al ccie c(q). 

1. Suppose that sj < g"'"(C*). Then by definition, x{S{Aj),a'^) = X]c'ec(C') ^("^(^fc)' "^)- ^° 



1(4. a') 



E 



x{T{i,a)\{{i,a)})+x{S{Al),a^) 



For each class CI G c(Cj), the minimum quantity '^a^y^i^x{T{i,a)\{{i,a)}) + x{S{Al.),a'^) is 

exactly Z{Cl,Sk,a^,a^)- As a result, for each class CI / C^,, let < Sk < q~^{Cl), and for class 
Cl,, let < Sk' < q^iCl) - 1: 



Z{Cl,Sk,aKa^). 



Thus, we can find out Z{Cj, sj, a^, a') by dynamic programming. 
2. Suppose that sj = q~^{Cj). Note that this time since the class will be "saturated", the term 
x{S{Aj),a^) does not get any positive values Xia, provided that a G C* n (lL!_^t ^ ^liat)- 
x{S{A'j),a^) = Ylciec{c^.)^(^(^k)^^^) ^^'^ ^^^^ implies that 



^ x(r(i,a)\{(i,a)})+x(5(yiy,a^ 



aGAl 



Let be the lowest ranking applicant that ranks higher than a^. Then for each class, C^ G c(Cj), 
the minimum quantity J2aeAi ^(^(^) '^)}) '^^i'^i^k)^^^) exactly Z{Cl,Sk,a^,a^)- Assuming 
that for each class Cl ^ Cl,, let < < q^{Cl), and let < s^' < q'^{Cl,) — 1, we have 



Z(Cj, Sj,a^,a^) = x{T{i, a^)) + mm^^.j-^^^^ ._i ^ Z{Cl, Sk, a^, i 
As before, this can be calculated by dynamic programming. 



□ 



Now choose the smallest Z{CpQ{i) — l,a^,a'l") among all possible who rank higher than 
and assume that is the corresponding tuple. It is easy to see that among all feasible tuples A^ of 
length Q{i) whose lowest ranking applicant is a^, the one has the smallest comb value x{K{i, S{A^)), 
is exactly the tuple AJ U {a^}. 

Separating Constraint ([5]) We again make use of dynamic programming. The idea is similar to 
the previous one and the task is much simpler, so we will be brief. 

Suppose that we are checking all the class-tuples ¥*• corresponding to class Cj. Let T*^^ C T*- be 

the subset of class-tuples whose lowest ranking applicants is . We need to find out the class-tuple 
t G T* t with the smallest value 
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r^(T(.,at)\{(^,at)}) + a)\{(., a)}), 

J, a' 

and check whether this value is no less than XlaGC'nL* -^ia- -'^^ then we are sure that all class- 

j -<a't 

tuples in T*^^ satisfy Constraint ([5]), otherwise, we find a violated constraint. The above quantity can 
be easily calculated by dynamic programming as before. 

E A Counter Example for Section 14.21 

The example shown in Figure U] contains five stable matchings. If we apply the median choice oper- 
ation on all of them, we get the stable matching /_i2, which does not give institutes ii and i2 their 
lexicographical median outcome. 



Institute Preferences 
i2:a^awa2aia4,as 



Classifications 

Cl = {ai,a2}, 6*2 = {as, 04} 
Cl = {ai, 02}, C| = {aa, 04} 



Class Bounds 

Q(il) = 2, q+{Cl) = l,q+{Cl) = l 
Q(i2) = 2, q+{Cl) = l,q+{Cl) = l 
Qih) = 4 



Applicant Preferences 

01:121113 

a2.iii2i3 

a3:«2ii«3 

a4:iii2i3 

ay-.isii 
02:4312 
0^^:1312 



Stable Matchings 

Ml = la;, fly), (12; flz, a™), (13; ai, 02, 03, 04)} 
M2 = {(*i;ai,a3), («2;a2,a4), («3; a^:, fly, Oz, a™)} 
M3 = {(n; ai,a4), (12; 02, as), (13; fla:, Oy, a2,a„)} 
/i4 = {(ii; 02,03), (12; oi, 04), (13; ffli, flj/, Oz,0„)} 

M5 = {(*i; 012,04), (22; 01, 03), (13; O3:, Oy, Oz,0„)} 

Fig. 4. An example of median choice stable matching which does not give the institutes their lexico- 
graphically median outcome. 



35 



